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A versatile method for combining density functional theory (DFT) in the local density approx- 
imation (LDA) with dynamical mean-field theory (DMFT) is presented. Starting from a general 
basis-independent formulation, we use Wannier functions as an interface between the two theories. 
These functions are used for the physical purpose of identifying the correlated orbitals in a spe- 
cific material, and also for the more technical purpose of interfacing DMFT with different kinds of 
band-structure methods (with three different techniques being used in the present work). We ex- 
plore and compare two distinct Wannier schemes, namely the maximally-localized- Wannier-function 
(MLWF) and the A'-th order muffin-tin-orbital (NMTO) methods. Two correlated materials with 
different degrees of structural and electronic complexity, SrVOs and BaVSa, are investigated as case 
studies. SrVOs belongs to the canonical class of correlated transition-metal oxides, and is chosen 
here as a test case in view of its simple structure and physical properties. In contrast, the sulfide 
BaVSs is known for its rich and complex physics, associated with strong correlation effects and 
low-dimensional characteristics. New insights into the physics associated with the metal-insulator 
transition of this compound are provided, particularly regarding correlation-induced modifications 
of its Fermi surface. Additionally, the necessary formalism for implementing self-consistency over 
the electronic charge density in a Wannier basis is discussed. 

PACS numbers: 71.30.+h, 71.15.Mb, 71.10.Fd, 75.30.Cr 



I. INTRODUCTION 

One of the fundamental points that underUes the rich 
physics of strongly correlated electron systems is the 
competition between the electrons tendency to localize, 
and their tendency to dclocalize by forming quasipar- 
ticle (QP) bands. Traditional effective single-particle 
(i.e. band-structme) theories emphasize the latter as- 
pect, which is appropriate when the kinetic energy dom- 
inates. For such materials, computational techniques 
based on electronic density functional theory (DFT) (see 
e.g. Refs. for reviews) have nowadays reached a very 
high degree of accuracy and yield remarkable agreement 
with experiment. 

In correlated materials however, the screened Coulomb 
interaction is a major aspect of the problem, which can- 
not be treated perturbatively, and independent-particle 
descriptions fail. Albeit the representability of the elec- 
tronic charge density by a set of Kohn-ShaniS. (KS) or- 
bitals is still guaranteed in most cases, this raises the 
question of whether such a representation is physically 
appropriate. Furthermore, the description of excited 
states of the many-particle system must be based on 
other observables than just the charge density, such as 
the energy-dependent spectral function. Any appropriate 
theoretical framework must then treat band formation 
(best described in momentum space) and the tendency 
to localization (best described in real space) on an equal 
footing. For this reason, there is an increasing aware- 
ness that many-body descriptions must also include real- 



space, orbitally resolved, descriptions of the solid, close 
to the quantum chemistry of the material under consid- 
eration^. In correlated metals, the coexistence of coher- 
ent QP bands at low energy with high-energy incoherent 
Hubbard bands (which originate from atomic-like transi- 
tions persisting in the solid state) is a vivid demonstra- 
tion that a dual description (both in momentum space 
and in real space) is needed. Such a dual description 
is at the heart of dyn amical mean-field theory (DMFT) 
(see e.g. Refs. |y|38||9||lO||ll| for reviews), which in re- 
cent years has proven to be a tool of choice for treating 
strong-correlation effects. This theory has been success- 
fully combined with electronic-structure methods within 
the framework of the local density approximation^ (LDA) 
to DFTiSJ^ (also labeled as LDA-hDMFT), or so-called 
GW formalismsili^. 

A central physical issue in merging the momentum- 
space and local descriptions within those many-body ap- 
proaches is the identification of a subspace of orbitals 
in which correlations are treated using non-perturbative 
many-body techniques. Furthermore, an important tech- 
nical issue is the choice of a convenient basis set for in- 
terfacing the many-body and the band-structure parts of 
the calculation. Because the original Wannier construc- 
tioriii is based on a decomposition of the extended Bloch 
states into a superposition of rather localized orbitals, it 
appears that appropriate generalizations of this construc- 
tion leading to well-localized basis functions should pro- 
vide an appropriate framework for many-body calcula- 
tions. Exploring this in detail is the main purpose of this 
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paper. In fact, there has been recently a growing activity 
associated with the Wannier formaUsm in the context of 
many-body approaches. The use of Wannier basis sets 
in the LDA+DMFT context has been lately pioneered 
by several groups, using either the A'^-th order muffin- 
tin orbital (NMTO) framcworh^?i^?i^9i^^ or other types 
of Wannier constructions based on the linear muffin-tin 
orbital (LMTO) £ramewory2i2i212^. For a detailed pre- 
sentation of such implementations see Ref. . Further- 
more, the computation of many-body interaction param- 
eters has also been discussed^SiSi. 

In this context, the main motivations of the present 
article arc the following: 

i) We give a presentation of the LDA-I-DMFT formal- 
ism in a way which should make it easier to inter- 
face it with a band-structure method of choice. To 
this aim, we are careful to distinguish between two 
key concepts: the orbitals defining the correlated 
subspacc in which a many-body treatment is done, 
and the specific basis set which is used in order 
to interface the calculation with a specific band- 
structure method. The LDA-I-DMFT approach is 
first presented in a manner which makes no ref- 
erence to a specific basis set, and then only some 
technical issues associated with choosing the basis 
set for implementation are discussed. 

ii) It is explained how the Wannier-functions for- 
malism provides an elegant solution to both the 
physical problem of identifying the correlated or- 
bitals and to the more technical issue of interfacing 
DMFT with basically any kind of band-structure 
method. So far the LDA-I-DMFT technique has 
been implemented with band-structure codes based 
on muffin-tin-orbital (MTO)-like representations^^. 
Although this realization is very successful, we feel 
that broadening the range of band-structure meth- 
ods that can be used in the LDA-I-DMFT context 
may make this method accessible to a larger part 
of the band-structure community, hence triggering 
further progress on a larger scale. As an example, 
one could think of problems involving local struc- 
tural relaxations, which are more difficult to handle 
within the MTO formalism than in plane-wave like 
approaches. 

iii) In this work, two different Wannier construc- 
tions are applied and the corresponding results 
are compared in detail. Though there are nu- 
merous ways of constructing Wannier(-like) func- 
tions we have chosen such methods that derive 
such functions in a post-processing step from a 
DFT calculation. In this way the method is, 
at least in principle, independent of the underly- 
ing band-structure code and therefore widely ac- 
cessible. First, we used the maximally-localized- 
Wannicr-functions (MLWFs) method proposed by 
Marzari, Vanderbilt and Souz a^^i'^" . Second, we 



constructed Wannier functions using the A^-th or- 
der MTO (NMTO) framework following Andersen 
and coworkcrsSiiSi^^ which has first been used in 
the LDA-hDMFT context in Re f. Ill and actively 
used since then (e.g. Refs. [lllllalli). Note that 
the NMTO method also works in principle with any 
given, not necessarily MTO-determined, KS effec- 
tive potential. However, in practice, this construc- 
tion is presently only available in an MTO environ- 
ment. 

iv) We also consider the issue of fully self-consistent 
calculations in which many-body effects are taken 
into account in the computation of the electronic 
charge density. Appendix ^ is devoted to a 
technical discussion of implementing charge self- 
consistency, with special attention to the use of 
Wannier basis sets also in this context. How- 
ever, the practical implementation of charge self- 
consistency in non-MTO based codes is ongoing 
work, to be discussed in detail in a future publi- 
cation. 

Two materials with correlated 3d electrons serve as 
testing grounds for the methods developed in this pa- 
per, namely the transition-metal oxide SrVOs and the 
sulfide BaVSs. Nominally, both compounds belong to 
the class of 3d^ systems, where due to crystal-field split- 
ting the single d electron is expected to occupy the t2g 
states only. The latter form partially filled bands in an 
LDA description. The two compounds have very differ- 
ent physics and exhibit different degrees of complexity in 
their electronic structure. The metallic pcrovskite SrVOs 
has perfect cubic symmetry over the temperature regime 
of interest and displays isolated i2g-like bands at the 
Fermi level, well-separated from bands higher and lower 
in energy. Its physical properties suggest that it is in 
a regime of intermediate strength of correlations. Many 
experimental results arc available for this material (for a 
detailed list see Sec. IIII A l|l and it has also been thor- 
oughly investigated theoretically in the LDA+DMFT 
frameworkiSii9,22,34,37,38,39^ Poj. ^j^ggg reasons, SrVOg 

is an ideal test case for methodological developments. 

In contrast, BaVSs is much more complex in both its 
electronic structure and physical properties. The sul- 
fide displays several second-order transitions with de- 
creasing temperature, including a metal-insulator transi- 
tion (MIT) at -70 K. Additionally, the low-energy LDA 
bands with strong t2g orbital character are entangled 
with other bands, mainly of dominant sulfur character, 
which renders a Wannier construction more challenging. 
In this paper, the Wannier-based formalism is used for 
BaVSa to investigate correlation-induced changes in or- 
bital populations, and most notably, correlation-induced 
changes in the shape of the different Fermi-surface sheets 
in the metallic regime above the MIT. In the end, these 
changes are key to a satisfactory description of the MIT. 

This article is organized as follows. Section ^ in- 
troduces the general theoretical formalism. First, the 
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LDA+DMFT approach is briefly reviewed in a way which 
does not emphasize a specific basis set. Then, the issue 
of choosing a basis set for implementation and interfac- 
ing DMFT with a specific band-structure method is dis- 
cussed. FinaUy, the Wannier construction is shown to 
provide an elegant solution for both picking the corre- 
lated orbitals and practical implementation. The differ- 
ent Wannier constructions used in this paper are briefly 
described, followed by some remarks on the calculational 
schemes employed in this work. In Sect. Illll the results 
for SrVOs and BaVSa are presented. To this aim we dis- 
cuss separately the LDA band structure, the correspond- 
ing Wannier basis sets and the respective LDA+DMFT 
results. Appendices are devoted to the basic formalism 
required to implement self-consistency over the charge 
density and total energy calculations, as well as further 
technical details on the DFT calculations. 



II. THEORETICAL FRAMEWORK 

A. Dynamical mean-field theory and electronic 
structure 

1. Projection onto localized orbitals 

Dynamical mean-field theory provides a general frame- 
work for electronic structure calculations of strongly cor- 
related materials. A main concept in this approach is a 
projection onto a set of spatially localized single-particle 
orbitals {|xRm)}, where the vector R labels a site in 
the (generally multi-atom) unit cell and m denotes the 
orbital degree of freedom. These orbitals generate a sub- 
space of the total Hilbert space, in which many-body 
effects will be treated in a non-perturbative manner. In 
the following, we shall therefore refer to this subspace as 
the "correlated subspace" C, and often make use of the 
projection operator onto this correlated subspace, defined 
as: 

p(f)^^|XR™)(xR™|. (1) 

For simplicity, we restrict the present discussion to the 
basic version of DMFT in which only a single correlated 
site is included in this projection. In cluster generaliza- 
tions of DMFT, a group of sites is taken into account. 
Also, it may be envisioned to generalize the method in 
such a way that R could stand for other physically des- 
ignated real-space entities (e.g. a bond, etc.). 

Because the many-body problem is considered in this 
projected subspace only, and because it is solved there in 
an approximate (though non-perturbative) manner, dif- 
ferent choices of these orbitals will in general lead to dif- 
ferent results. How to properly choose these orbitals is 
therefore a key question. Ultimately, one might consider 
a variational principle which dictates an optimal choice 
(cf. Appendix IB|) . At the present stage however, the 
guiding principles are usually physical intuition based on 



the quantum chemistry of the investigated material, as 
well as practical considerations. Many early implemen- 
tations of the LDA-I-DMFT approach have used a linear 
muffin-tin orbital^^jXMTO) basis for the correlated or- 
bitals (e.g. Refs. [l3i.4Qi |'l. This is very natural since in 
this framework it is easy to select the correlated subspace 
C regarding the orbital character of the basis functions: 
e.g., d character in a transition-metal oxide, / charac- 
ter in rare-earth materials, etc.. The index m then runs 
over the symmetry-adapted basis functions (or possibly 
the "heads" of these LMTOs) corresponding to this se- 
lected orbital character. Exploring other choices based 
on different Wannier constructions is the purpose of the 
present paper. In this context, the index m should be 
understood as a mere label of the orbitals spanning the 
correlated subset. For simplicity, we shall assume in the 
following that the correlated orbitals form an orthonor- 
mal set: (xR„|xR'm')=<5RR'^mm'- This may not be an 
optimal choice for the DMFT approximation however, 
which is better when interactions are more local. Gen- 
eralization to non-orthogonal sets is yet straightforward 
by introducing an overlap matrix (see e.g. Ref. [^). 

2. Local observables 

There are two central observables in the LDA+DMFT 
approach to electronic structure. The first, as in DFT, 
is the total electronic charge density p{r). The second 
is the local one-particle Green's function Gioc(*w„) pro- 
jected onto C, with components GRm.Rm' (*i^n)- Both 
quantities are related to the full Green's function of the 
solid G(r, r';ia;„) by: 

p(r)^i^G(r,r;»o.„)e--°^ (2) 

n 

G'Zn'i^^n) = I ^drdr'xLn(r-R)XR™'(r'-R)G(r,rVo.„)- 

The last expression can be abbreviated as a projection of 
the full Green's function operator G according to 

Gioc^P^^GP^K (3) 

In these expressions, we have used (for convenience) the 
Matsubara finite-temperature formalism, with ujn={2n + 
1)tt/P and /3=1/A:bT. The Matsubara frequencies are re- 
lated via Fourier transformation to the imaginary times 
T. Note that the factor e*""° in (01 ensures the con- 
vergence of the Matsubara sum which otherwise falls of 
as l/ujn- We have assumed, for simplicity, that there is 
only one inequivalent correlated atom in the unit cell, so 
that Gioc does not carry an atom index (generalization is 
straightforward) . In the following we will drop the index 
R if not explicitly needed. 

Taking the KS Green's function Gks as a reference, 
the full Green's function of the solid can be written in 
operator form, as: G^^=Gj^g — AS, or more explicitly 
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(atomic units are used throughout this paper): 



G(r, r';iw„) 



lUJn+ Vks - AS 



(4) 

Here n is the chemical potential and Vks the KS effective 
potential, which reads: 



|r-N,| J |r-r'| Spir) 
= 14xt(r) + V%(r) + 14c(r), (5) 

where Q, N are the charges, lattice vectors of the atomic 
nuclei, i runs over the lattice sites, V^xt is the exter- 
nal potential due to the nuclei, Vr denotes the Hartree 
potential and Vxc is the exchange-correlation potential, 
obtained from a functional derivative of the exchange- 
correlation energy E^c- For the latter, the LDA (or 
generalized-gradient approximations) may be used. Re- 
call that Vks (r) is determined by the true self-consistent 
electronic charge density given by Eq. (0) (which is mod- 
ified by correlation effects, and hence differs in general 
from its LDA value, see below). 

The operatorAE in Eq. describes a many-body 
(frequency-dependent) self-energy correction. In the 
DMFT approach, this self-energy correction is con- 
structed in two steps. First, AS is derived from an ef- 
fective local problemii (or "effective quantum impurity 
model" ) within the correlated subspace C via: 



AS, 



v^imp 



,{iU!n) 



(6) 



whereby Sdc is a double-counting term that corrects for 
correlation effects already included in conventional DFT. 
The self-energy correction to be used in and subse- 
quently in (|12I13|I . is then obtained by promoting © to 
the lattice, i.e., 

AS(r,r';itj„) = 
J2 X™(r-R-T)x„,(r'-R-T)AS„„„,(za;„), (7) 

Tmm' 

where T denotes a direct lattice vector. The key ap- 
proximation is that the self-energy correction is non-zero 
only inside the (lattice-translated) correlated subspace, 
i.e., AE=AE(''', hence exhibits only on-site components 
in the chosen orbital set. 



3. Effective quantum impurity problem 

The local impurity problem can be viewed as an ef- 
fective atom involving the correlated orbitals, coupled to 
a self-consistent energy-dependent bath. It can be for- 
mulated either in Hamiltonian form, by explicitly intro- 
ducing the degrees of freedom of the effective bath, or as 
an effective action in which the bath degrees of freedom 



have been integrated out. In the latter formulation, the 
action of the effective impurity model reads: 



5'imp = - / / dTdr'Y^ 4i^(T)[^o \un'{T - T')d„^,„{T') 



+ I drHij [{d\^^\ dyna}) 



(8) 



In this expression, d^^ are the Grassmann vari- 

ables corresponding to orbital Xm for spin cr, Hu is a 
many-body interaction (to be discussed in scction lTl G 2|l . 
and QQ[iijJn) is the dynamical mean-field, determined self- 
consistently (see below), which encodes the coupling of 
the embedded atom to the effective bath. This quantity is 
the natural generalization to quantum many-body prob- 
lems of the Weiss mean-field of classical statistical me- 
chanics. Its frequency dependence is the essential feature 
which renders DMFT distinct from static approaches 
such as e.g. the LDA+U method^. The frequency de- 
pendence allows for the inlcusion of all (local) quantum 
fluctuations. Thereby the relevant (possibly multiple) 
energy scales are properly taken into account, as well as 
the description of the transfers of spectral weight. One 
should note that the dynamical mean- field Qo{iijJn) for- 
mally appears as the bare propagator in the definition 
of the effective action for the impurity ©. However, 
its actual value is only known after iteration of a self- 
consistency cycle (detailed below) and hence depends on 
many-body effects for the material under consideration. 

The self-energy correction is obtained from the impu- 
rity model as: 



Simp(iw„) = 5q ^(iw„) - GiJp(iw„) 



(9) 



in which Gi,„p is the impurity model Green's function, 
associated with the effective action ^ and defined as: 



^imp / 

mm' V 



')^-(td_(r)dl,,,(r'))s. 



(10) 



where T stands for time-ordering. Note that computing 
this Green's function, given a specific Weiss dynamical 
mean- field GoiioJn) is in fact the most demanding step in 
the solution of the DMFT equations. 



4- Self- consistency conditions 

In order to have a full set of self-consistent equations, 
one still needs to relate the effective impurity problem 
to the whole solid. Obviously, the dynamical mean-field 
Qo{iu!n) is the relevant link, but we have not yet specified 
how to determine it. The central point in DMFT is to 
evaluate Qo{iuJn) in a self-consistent manner by request- 
ing that the impurity Green's function coincides with the 
components of the lattice Green's function projected onto 
the correlated subspace C, namely that: 



(11) 
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(Appendix A) 
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FIG. 1; Complete self-consistency loop for LDA-fDMFT. The charge density p determines the KS potential Vks, from which 
KS eigenvalues ekf and eigenfunctions i/j^^ follow. The KS Green's function is then constructed and passed on to the DMFT 
cycle (in practice, the KS Hamiltonian Hks is constructed in the basis set used to implement the method, and transferred to the 
DMFT cycle). The DMFT loop consists in i) solving the effective impurity problem for the impurity Green's function, hence 
obtaining an impurity self-energy, ii) combining the self-energy correction with the KS Green's function in order to obtain the 
local Green's function Gioc projected in the correlated subset and iii) obtaining an updated Weiss mean-field. An initial guess 
for the Weiss dynamical mean-field must be made at the beginning of the DMFT loop, e.g. by choosing So"'* = P^'^'^GksP^'^'' ■ 
Is the DMFT loop converged, the chemical potential is updated in order to ensure global charge neutrality, and the new 
charge density (including many-body effects) is constructed (described in Appendix. EJ- This new density determines a new 
KS potential. Note that in addition one may want to update the set {|xm)} when preparing for the next DMFT loop (cf. 
AppendixIQ. The whole process must be iterated until the charge density, the impurity self-energy and the chemical potential 
are converged. In practice, good convergence of the DMFT loop is reached before a new p is calculated. Note that in the 
present paper using a Wannier implementation, the global self-consistency on the charge density is not implemented in practice. 
Thus the self-consistent LDA Hamiltonian Hks enters the DMFT loop, which is iterated until convergence of the self-energy. 



or, in explicit form, using H2I4I6|) and (Q: 

GZZ'i^^n) ^ jj drdr'x:,(r-R)x„,,(r'-R) x (12) 



X ( r 



2 



iuJn + /i + Fks — ( Sjmp — Sdc 



In this representation, it is clear that the self-consistency 
condition involves only impurity quantities, and therefore 
yields a relation between the dynamical mean-field 
and Gimp which, together with the solution of the impu- 
rity problem Eqs. (|8I10|) fully determines both quantities 
in a self-consistent way. 

The effective impurity problem Q can in fact be 
thought of as a reference system allowing one to repre- 
sent the local Green's function. This notion of a reference 
system is analogous to the KS construction, in which the 
charge density is represented as the solution of a single- 



electron problem in an effective potential (with the dif- 
ference that here, the reference system is an interacting 
one). 

Finally, by combining Q and (0J the electronic charge 
density is related to the KS potential by: 




V2 ^ 

lujn + /-i + Vks 



AS 



(13) 

Expression (|13|) calls for two remarks. Firstly, many- 
body effects in C affect via AS the determination 
of the charge density, which will thus differ at self- 
consistency from its LDA value. Secondly, the fa- 
miliar KS representation of p(r) in terms of virtu- 
ally independent electrons in an effective static poten- 
tial is modified in LDA+DMFT in favor of a non- 
local and energy-dependent (retarded) potential given by 
Vks(r)^(r- r')(5(T-r')+AS(r,r';T-T'). In Appendix|l| 
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we give a more detailed discussion of the technical aspects 
involved in calculating the charge density from expres- 
sion 1)13(1 . However, we have not yet implemented this 
calculation in practice in our Wannier-based code: the 
computations presented in this paper are performed for 
the converged p(r) obtained at the LDA level. Finally, let 
us mention that the LDA+DMFT formalism and equa- 
tions presented above can be derived from a (free-energy) 
functional^ of both the charge density and the projected 
local Green's function, T[p, Giod- This is reviewed in Ap- 
pendix where the corresponding formula for the total 
energy is also discussed. 

Fig. n gives a synthetic overview of the key steps in- 
volved in performing a fully self-consistent LDA-I-DMFT 
calculation, irrespective of the specific basis set and band- 
structure code chosen to implement the method. 



5. Double-counting correction 

We briefly want to comment on the double-counting 
(DC) correction term. Since electronic correlations are 
already partially taken into account within the DFT ap- 
proach through the LDA/GGA exchange-correlation po- 
tential, the double-counting correction Sdc has to cor- 
rect for this in LDA-I-DMFT. The problem of precisely 
defining DC is hard to solve in the framework of con- 
ventional DFTil^. Indeed, DFT is not an orbitally- 
rcsolved theory and furthermore the LDA/GGA does not 
have a diagrammatic interpretation (like simple Hartree- 
Fock) which would allow to subtract the corresponding 
terms from the DMFT many-body correction. Simply 
substracting the matrix elements of Vh and V^c in the 
correlated orbital subset C from the KS Green's function 
to which the many-body self-energy is applied to is not 
a physically reasonable strategy. Indeed, the DMFT ap- 
proach (with a static, frequency-independent Hubbard 
interaction) is meant to treat the low-energy, screened 
interaction, so that the Hartree approximation is not an 
appropriate starting point. Instead, one wants to bene- 
fit from the spatially-resolved screening effects which are 
already partially captured at the LDA level. In prac- 
tice, the DC terms introduced for LDA-I-U, i.e., "fully- 
localized limit"— and "around mean Reld"^^^, appear 
to be reasonable also in the LDA-I-DMFT framework. It 
was recently shown^^, that the fuUy-localized-limit form 
can be derived from the demand for discontinuity of the 
DFT exchange-correlation potential at integer filling. 

The DC issue in fact has a better chance to be re- 
solved in a satisfactory manner, from both the physical 
and formal points of view, when the concept of local in- 
teraction parameters is extended to frequency-dependent 
quantities (e.g. a frequency-dependent Hubbard interac- 
tion U{ll!)), varying from the bare unscreened value at 
high frequency to a screened value at low energy, and 
determined from first principles. The GW+DMFT con- 
struction, and the extended DMFT framework, in which 
this quantity plays a central role and is determined self- 



consistently on the same footing as the one-particle dy- 
namical mean field, may prove to be a fruitful approach 
in this respect. 



6. Implementation: choice of basis sets and Hilbert spaces 

In the previous sections, care has been taken to write 
the basic equations of the LDA+DMFT formalism in a 
basis-independent manner. In this section, we express 
these equations in a general basis set, which is essen- 
tial for practical implementations, and discuss advan- 
tages and drawbacks of different choices for the basis 
set. At this point, a word of caution is in order: it is 
important to clearly distinguish between the set of local 
orbitals {|Xm)} which specifies the correlated subspace, 
and the basis Junctions which one will have to use in order 
to implement the method in practice within a electronic- 
structure code. Different choices for {|xm)} will lead to 
different results, since DMFT involves a local approxima- 
tion which has a different degree of accuracy in diverse 
orbital sets. In contrast, once the correlated orbital set 
is fixed, any choice of basis set can be used in order to 
implement the method, with in principle identical results. 

Let us denote a general basis set by {|i3ka)}, in which 
k runs over the Brillouin zone (BZ) and a is a label for 
the basis functions. For example if the KS (Bloch) wave 
functions are used as a basis set, a=v is a band index and 
the basis functions are '0ki/(r)=e'''''"uki/(r)- In the case 
of a pure plane- wave basis set, a~G runs over reciprocal 
lattice vectors. For a Bloch sum of LMTOs q;=(R^to) 
runs over sites in the primitive cell and orbital (angular 
momentum) quantum number. For hybrid basis sets one 
may have a=(RZm, G). As an example for the latter 
serves the linear augmcnted-plane-wave^^ (LAPW) basis 
set, though here in the end it is summed over the orbital 
indices and the basis is finally labelled by G only. 

Consider now the DMFT self-consistency condition 
In the (yet arbitrary) basis set {l-Bka)}; its explicit 
expression in reciprocal space reads (correct normaliza- 
tion of the fc-sum is understood): 

- EE(x™i^k„)(i?ka'ix!;')x 



{I. 



H 



KS 



(k)- AS(k,io.„)]-^) (14) 



where |Xm)= X)t e''' '^'^+^^ \Xm) denotes the Bloch trans- 
form of the local orbitals. In this expression, HKs(k) is 
the KS Hamiltonian at a given fc-point, expressed in the 
basis set of interest: 



7?Ks(k) -El^ka)(Ska' 



E ^ki. (-BkQ I ) ( V'ki/ 1 Ska' 

(15)" 

with {skiy, V'kiy} the set of KS eigenvalues and wave func- 
tions: 

V2 - 1 

-^ + Vks iV'k^) = eki.l^/'ki.) (16) 
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The self-energy correction, in the chosen basis set, reads: 



(17) 



and it should be noted that, although purely local when 
expressed in the set of correlated orbitals, it acquires in 
general momentum-dependence when expressed in an ar- 
bitrary basis set. 

The self-consistency condition (|14|) is a central step in 
interfacing DMFT with a chosen band-structure method. 
Given a charge density p(r), the effective potential Vks (r) 
is constructed, and the corresponding KS equations Hlt)|) 
are solved (Fig. in a manner which depends on the 
band-structure method. Each specific technique makes 
use of a specific basis set {|i3ka)}- The KS Hamiltonian 
serves as an input to the DMFT calculation for C, which 
is used in (|14|l to recalculate a new local Green's function 
from the impurity self-energy, and hence a new dynamical 
mean-field from P(7^=G[~J -I- Simp- 

A remark which is important for practical implemen- 
tation must now be made. Although Gioc(i'j-'„), i.e., the 
right-hand side of (|14|l . can be evaluated in principle 
within any basis set {|i?ka)}, the computational effort 
may vary dramatically depending on the number Nb of 
basis functions in the set. According to H14|l . this com- 
putation involves an inversion of a Nb x Nb matrix at 
each k-point and at each frequency iuJm followed by a 
summation over fc-points for each frequency. Since the 
number of discrete frequencies is usually of the order of 
a few thousands, this procedure is surely feasible within 
a minimal basis set such as, e.g., LMTOs. In the latter 
case, the correlated orbitals may furthermore be chosen 
as a specific subset of the basis functions (e.g. with d 
character in a transition-metal oxide) -or possibly as the 
normalized "heads" corresponding to this subset-, mak- 
ing such basis sets quite naturally tailored to the prob- 
lem. In contrast, computational efficiency is harder to 
reach for plane-wave like basis sets in the LDA-I-DMFT 
context. For such large basis sets, the frequency depen- 
dence substantially increases the already large numerical 
effort involved in static schemes such as LDA or LDA-I-U. 
Furthermore, another more physical issue in using plane- 
wave based codes in the DMFT context is how to choose 
the local orbitals {|Xm)} which define the correlated sub- 
set. Because the free-electron like basis functions usually 
do not have a direct physical connection to the quantum 
chemistry of the material problem at hand, these orbitals 
must be chosen quite independently from the basis set it- 
self. 

To summarize, when implementing LDA+DMFT in 
practice, a decision must be made on the following two 
issues: 

-i) The first issue is a physical one, namely how to 
choose the local orbitals Xm spanning the correlated sub- 
space C. The quality of the DMFT approximation will in 
general depend on the choice of C, and different choices 



may lead to different results. Obviously, one would like 
to pick C in such a way that the DMFT approximation is 
better justified, which is intuitively associated with well- 
localized orbitals. 

-ii) The second point is a technical, albeit important, 
one. It is the choice of basis functions {|i?kQ}} used for 
implementing the self-consistency condition H14|) . As dis- 
cussed above and as clear from (|14|) . computational effi- 
ciency requires that as many matrix elements (-BkalXm) 
as possible arc zero (or very small), i.e., such that Xm 
has overlap with only few basis functions. 

As discussed above, both issues demand particular at- 
tention when using band-structure methods based on 
plane-wave techniques, because those methods do not 
come with an obvious choice for the orbitals Xm ^'^'^ be- 
cause the demand for well-localized Xm implies that they 
will overlap with a very large number of plane waves. 

In this paper, we explore the use of WFs as an elegant 
way of addressing both issues i) and ii), leading to a 
convenient and efficient interfacing of DMFT with any 
kind of band-structure method. 



B. Wannier functions and DMFT 

1. General framework and Wannier basics 

Let us outline the general strategy that may be used 
for implementing LDA+DMFT using Wannier functions 
(WFs), postponing technical details to later in this sec- 
tion. First, it is important to realize that a Wannier con- 
struction needs not be applied to all Bloch bands span- 
ning the full Hilbert space, but only to a smaller set W 
corresponding to a certain energy range, defining a subset 
of valence bands relevant to the material under consid- 
eration. To be concrete, in a transition-metal oxide for 
example, it may be advisable to keep bands with oxygen 
2p and transition-metal Sd character in the valence set 
yy. WFs spanning the set W may be obtained by per- 
forming a (k-dcpendent) unitary transformation on the 
selected set of Bloch functions. This unitary transforma- 
tion should ensure a strongly localized character of the 
emerging WFs. Among the localized WFs spanning W, 
a subset is selected which defines the correlated subspace 
CcyV. For transition-metal oxides, C will in general cor- 
respond to the WFs with d character. 

The correlated orbitals Xm thus identified with a 
certain set of WFs generating C. It is then recommend- 
able (albeit not compulsory) to choose the (in general 
larger) set of WFs generating the valence set W, as ba- 
sis functions in which to implement the self-consistency 
condition H14|) . Indeed, the KS Hamiltonian can then be 
written as a matrix with diagonal entries corresponding 
to Bloch bands outside W, and only one non-diagonal 
block corresponding to W. It follows that the self- 
consistency condition H14|l may be expressed in a form 
which involves only the knowledge of the KS Hamiltonian 
within W and requires only a matrix inversion within this 
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subspace, as detailed below. Hence, using WFs is an ele- 
gant answer to both points i) and ii) above: it allows to 
build correlated orbitals defining the set C with tailored 
localization properties, and by construction only the ma- 
trix elements (Xmlw^a) with Wa a WF in the set W arc 
non-zero. 

We now describe in more details how WFs are con- 
structed. Within the Born-von Karman periodic bound- 
ary conditions, the effective single-particle description 
of the electronic structure is usually based on extended 
Blocli functions '(/'ki/i which are classified with two quan- 
tum numbers, the band index v and the crystal momen- 
tum k. An alternative description can be derived in terms 
of localized WFsii, which are defined in real space via 
an unitary transformation performed on the Bloch func- 
tions. They arc also labeled with two quantum numbers: 
the index a which describes orbital character and po- 
sition, as well as the direct lattice vector T, indicating 
the unit cell they belong to. The relation between WFs 
and Bloch functions can be considered as the generaliza- 
tion to solids of the relation between "Boys orbitals"— 
and localized molecular orbitals for finite systems. It is 
crucial to realize, that the unitary transformation is not 
unique. In the case of an isolated band in one dimen- 
sion, this was emphasized long ago by W. Kohn^. He 
stated that infinitely many WFs can be constructed by 
introducing a k-depcndcnt phase </5(k), yet there is only 
one real high-symmetry WF that falls off exponentially. 
Hence in general </9(k) may be optimized in order to im- 
prove the spatial localization of the WF in realistic cases. 
This observation was generalized and put in practice for 
a group of several bands in Ref . . 

Let us consider the previously defined group W of 
bands of interest. A general set of WFs corresponding 
to this group can be constructed as^S 



t«a(r-T) 



V 



BZ 



dke-'^-T ^ U^^J^^,{v), (18) 



V denoting the volume of the primitive cell. The WF 
(r|u>Ta) only depends on r— T, since ■i/)jj.^,(r)=c*'' ''Mki/(r); 
with Ukiy a periodic function on the lattice. The uni- 
tary matrix reflects the fact that, in addition to 
the gauge freedom with respect to a k-dependent phase, 
there is the possibility of unitary mixing of several crys- 
tal wave functions in the determination of a desired WF. 
Optimization of these degrees of freedom allows one to 
enforce certain properties on the WFs, including the de- 
mand for maximal localization (see next paragraph). Of 
course, the extent of the WF still depends on the spe- 
cific material problem. Due to the orthonormality of the 
Bloch functions, the WFs also form an orthonormal basis: 
{w'Ta\wTa')— ^TT^aa'- Morc On the general properties 
and specific details of these functions ma y be foun d in 
the original literatureiLiLiSia, or Refs. |23l2ii3(H and 
references therein. 

Here, LDA-I-DMFT will be implemented by selecting 
a certain subset {wm, m S C} of the WFs {wa, a £ W} 



as generating the correlated subset. Thus we directly 
identify {|x,„)} with a specific set of WFs. Note again 
that this is a certain choice, and that other choices are 
possible (such as identifying {Ix^}} from only parts of 
the full WFs through a projection). With our choice, the 
functions 



T 



E 



C/^^.Vk.), (19) 



will be used in order to express the KS Hamiltonian and 
to implement the self-consistency condition H14|l . Be- 
cause the unitary transformation acts only inside W, 
only the block of the KS Hamiltonian corresponding to 
this subspace needs to be considered when implementing 
the self-consistency condition, hence leading to a quite 
economical and well-defined implementation. The KS 
Hamiltonian in the space W reads: 



(20) 



act' 6 W 

F„„-(k)= 
while the self-energy correction reads: 

AEt'^) = ^ [S""P(ic^„) - T.^%am' I^k™)(l«k™' I . 

7nm' G C k 

(21) 

Accordingly, the DMFT self-consistency condition takes 
the form: 



1, ^ 



(22) 



H|,T(k) 



In this expression, the matrix inversion has to be done for 
the full W-space matrix, while only the block correspond- 
ing to C has to be summed over k in order to produce the 
local Green's function Gioc in the correlated subspace, 
i.e. the r.h.s of (|22|) . In practice, the latter is inverted 
and added to Simp in order to produce an updated dy- 
namical mean-field according to: ^(^^^Gj^J + Sinip. This 
new dynamical mean-field is injected into the impurity 
solver, and the iteration of this DMFT loop leads to a 
converged solution of (|22|l (cf. Fig. 

In all the above, we have been careful to distinguish 
the (larger) space W in which the Wannier construction 
is performed, and the (smaller) subset C generated by the 
Wannier functions associated with correlated states. In 
some cases however, it may be possible to work within 
an energy window encompassing only the "correlated" 
bands (e.g when they are well separated from all other 
bands), and choose W=C. This of course leads to more 
extended Wannier functions than when the Wannier con- 
struction is made in a larger energy window. For the two 
materials considered in this paper, we shall nonetheless 
adopt this "massive downfolding" route, and work with 
yV=C. For the correlated perovskite SrVOa, the bands 
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originating from the ligand orbitals are well separated 
from the transition-metal ones. In other words, the size 
of the many-body interaction, say the Hubbard U, is ex- 
pected to be significantly smaller than the former level 
separation. In that case the minimal choice of a sub- 
set W=C involving only the d-likc WFs of the t2g panel 
is quite natural (see below). The situation is more in- 
volved for the BaVSa compound, since the S(3p) bands 
are strongly entangled with the t2g bands. Despite this 
stronger hybridization, it is not expected that the S(3p)- 
V(3d) level separation is the relevant energy scale, but 
still U. Hence we continue to concentrate on a disen- 
tangled t2g-likc panel, thereby integrating out explicit 
sulfur degrees of freedom. The resulting minimal basis 
is only "Wannier-like" , but nonetheless should provide a 
meaningful description of the low-energy sector of this 
material. It should be kept in mind however that the 
minimal choice W=C may become a critical approxima- 
tion at some point. For late transition- metal oxides in 
particular, the fact that p- and d-like bands are rather 
close in energy almost certainly implies that W must re- 
tain 0{2p) states, as well as transition-metal 3d states 
(while C will involve the 3d states only)^. 

2. Maximally-localized Wannier functions 

The maximally-localized Wannier function922i2£ (ML- 
WFs) are directly based on Eq. (|18|) . In order to ensure 
a maximally-localized Wannicr(-like) basis, the unitary 

fk) 

matrix Uav is obtained from a minimization of the sum 
of the quadratic spreads of the Wannier probability dis- 
tributions, defined as 

l]^^((r2)„-(r)i) , (0)„= fdrOMr)]'. 

a •' 

(23) 

Thus the quantity f2 may be understood as a functional 
of the Wannier basis set, i.e., 0=ri[{ii;Q}]. Starting from 
some inital guess for the Wannier basis, the formalism 
uses steepest-decent or conjugate-gradient methods to 
optimize C/i^'. Thereby, the gradient of Q. is expressed 
in reciprocal space with the help of the overlap matrix 

M^!^;'"' = K.i«k+q.') , (24) 

where q is connecting k vectors on a chosen mesh in recip- 
rocal space. Hence this scheme needs as an input the KS 
Bloch eigenfunctions V'ki^' '^'^ rather their periodical part 
Wki/. In the formalism, all relevant observables may be 
written in terms of m'^]'^'^ . The resulting MLWFs turn 
out to be real functions, although there is no available 
general proof for this property. 

In the following, two cases of interest shall be sepa- 
rately discussed: 

a. Bands of interest form a group of isolated bands. 
This is the case e.g for SrVOs discussed in this paper. 
The matrix M^''''^)^''^ has to be initially calculated from 



the KS Bloch eigenvectors V'kiyi where v runs over the 
bands defining the isolated group. Starting from 
according to the initial Wannier guess, the unitary trans- 
formation matrix will be updated iterativelySS. Corre- 
spondingly, the M matrices evolve as 

]V[(k,q) ^ u('')'fM^'''''^'^°^U^''+''^ . (25) 

The minimization procedure not only determines the in- 
dividual spreads of the WFs, but also their respective 
centers. Thus generally the centers do not have to coin- 
cide with the lattice sites as in most tight-binding rep- 
resentations. For instance, performing this Wannier con- 
struction for the four valence bands of silicon leads to 
WFs which are exactly centered in between the atoms 
along the bonding axesSS. 

h. Bands of interest are entangled with other hands. 
The handling of BaVSs discussed in this paper falls into 
this category. This case is not so straightforward, since 
before evaluating the MLWFs one has to decide on the 
specific bands subject to the Wannier construction. Lets 
assume there are N}, target bands, e.g. a t23-like mani- 
fold, strongly hybridized with iV^ other bands of mainly 
different character, e.g. s- or p-like bands. Then first the 
matrix (|24|l has to be calculated initially for the enlarged 
set of Nb+Nl^ bands. Within the latter set, the orbital 
character corresponding to the aimed at WFs may jump 
significantly. Thus new effective bands, associated with 
eigenvectors ip-^^, have to be constructed in the energy 
window of interest according to a physically meaningful 
description. 

To this aim, the functional ri[{wQ}] was decomposed 
in Ref. |3Cll| into two non- negative contributions, i.e., 
ri=r2/ -I- f2. Here ^2/ describes the spillage^ of the WFs 
between different regions in reciprocal space. The second 
part Cl measures to what extent the MLWFs fail to be 
eigenfunctions of the band-projected position operators. 
In the case of an isolated set of bands flj is gauge invari- 
ant. However it plays a major role in the case of entan- 
gled bands^, since here it may define a guiding quantity 
for "downfolding" the maximally (A^;,-f 7V^)-dimensional 
Hilbert space at each fc-point to a corresponding Hilbert 
space with maximal dimension A^^. The reason for this 
is that an initial minimization of Qj provides effective 
target bands with the property of "global smoothness of 
connection"—. Since flj measures the spillage, minimiz- 
ing it corresponds to choosing paths in reciprocal space 
with minimal mismatch within the reduced set of N^. In 
a second step Q is minimized for these effective bands, 
corresponding to the "traditional" procedure outlined for 
the isolated-bands case. Hence is now applied to the 
V'kiy- Note however that no true WFs in the sense of H18|) 
result from this procedure due to the intermediate cre- 
ation of effective bands. Yet the obtained Wannier-like 
functions are still orthonormal and stem from Bloch-like 
functions. 
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3. N-th order muffin-tin-orbital Wannier functions 

In this paper, we also consider another estabhshed 
route for the construction of locahzed Wannier(-hke) 
functions, namely the iV-th order muffin-tin-orbital 
(NMTO) method2ii^2i2a. This method is the latest 
development of the linear muffin-tin-orbital (LMTO) 
method^SiS. It uses multiple-scattering theory for an 
overlapping muffin-tin potential to construct a local- 
orbital minimal basis set, chosen to be exact at some 
mesh of iV-|-l energies, eo,...,eAr. This NMTO set is 
therefore a polynomial approximation (PA) in the energy 
variable to the Hilbert space formed by all solutions of 
Schrodinger's equation for an effective single-particle po- 
tential. In the present case this potential is given by the 
overlapping muffin-tin approximation to the KS potential 
(Eq. jSJ). Hence in contrast to the maximally-localized 
procedure, the NMTO-WFs for correlated bands may be 
generated without explicitly calculating the correspond- 
ing Bloch hmctions. 

Apart from its energy mesh, an NMTO set is specified 
by its members {R^m}, where Im denotes an angular- 
momentum character around site R, in any primitive cell. 
The RZm-NMTO is thus centered mainly at R and has 
mainly Im character. Moreover, for the NMTO set to be 
complete for the energies on the mesh, each NMTO must 
be constructed in such a way that its projections onto 
the RZm-channels not belonging to the {R^m}-set are 
regular solutions of Schrodinger's equationi24. Finally, 
in order to confine the R^m-NMTO, it is constructed in 
such a way that its projections onto all other channels 
belonging to the {RZm}-set vanish. 

For examploiS, the three isolated t2g bands of cubic 
SrVOa are spanned quite accurately by the quadratic 
(iV=2) muffin-tin orbital set which consists of the three 
(congruent) dxy, dyz and dxz NMTOs placed on each V 
site in the crystal. Locally, the dxy orbital has xy char- 
acter as well as minute other characters compatible with 
the local symmetry, but no yz or xz characters. On the 
O sites, the V dxy orbital has antibonding 0{p) and other 
characters compatible with the energy and the symme- 
try, in particular px character on O along the y-axis and 
Py character on the O along the x axis. On the Sr sites, 
there are small contributions which bond to 0(p). Fi- 
nally, on the other V sites, there can be no t2g character, 
but minute other characters are allowed by the local sym- 
metry. Note that when the symmetry is lowered, as is the 
case for the distorted pcrovskites CaVOs, LaTiOa, and 
YTiOs, there are less symmetry restrictions on the down- 
folded channels and the cation character of the V or Ti 
t2g NMTOs will increaseiSiiS. This describes a measur- 
able effect of cation covalency, and is not an artefact of 
the NMTO construction. 

The main steps in the NMTO construction arc thus: 
(a) numerical solution of the radial Schrodinger (or 
Dirac) equation for each energy on the mesh and for 
each I channel with a non-zero phase shift; (b) screening 
(or downfolding) transformation of the Korringa-Kohn- 



Rostocker (KKR) matrix for each energy on the mesh; 
and (c) formation of divided differences on the mesh of 
the inverse screened matrix to form the Lagrange matrix 
of the PA, as well as the Hamiltonian and overlap matri- 
ces in the NMTO representation. It should be noted that 
this procedure of downfolding plus PA differs from stan- 
dard Lowdin downfolding^ and is more accurate when 
iV>l. 

For an isolated set of bands and with an energy mesh 
spanning these bands, the NMTO set converges fast with 
N . The converged set spans the same Hilbert space as 
any Wannier set, and may even be more localized because 
the NMTO set is not forced to be orthonormal. Sym- 
metrical orthonormalization of the converged NMTO set 
yields a set of WFs w-Rim, which are atom-centered and 
localized. However this does not imply that the centre of 
gravity is the centre of the atom (see e.g. Fig. 5 and 6 of 
Ref. [iil). Note that NMTO-WFs have not been chosen 
to minimize the spread {wjiiml'r'^ IwRim) i but to satisfy 
the above-mentioned criterion of confinement. Using lo- 
calized NMTOs, it does not require a major computa- 
tional effort to form linear combinations which maximize 
any other suitable measure of localization. 



C. Calculational scheme 

1. Band-structure calculations and Wannier construction 

In the following we briefly name the different first- 
principles techniques that were used in the DFT part 
of the work. More technical details on the specific setups 
may be found in Appendix lUl 

The MLWF scheme was interfaced in this work with 
a mixed-basis^ pseudopotential^&S (MBPP) codeS. 
This band-structure program utilizes scalar-relativistic 
norm-conserving pseudopotentialaS and a basis of plane 
waves supplemented by non-overlapping localized func- 
tions centered at appropriate atomic sites. The localized 
functions, usually atomic functions for a given reference 
configuration, are necessary to allow for a reasonable 
plane-wave cutoff when treating electronic states with 
substantial local character. No shape approximations to 
the potential or the charge density are introduced and no 
MT spheres are utilized in this formalism. 

In addition, we also interfaced an already ex- 
isting MLWF schemaS£ with the all-electron, full- 
potential- linearized-augmented-plane- wave (FLAPW) 
method^SiSiiS^. This technique is fully self-consistent, 
i.e., all electrons are treated within the self-consistency 
procedure, and no shape approximations are made for 
the charge density and the potential. The core electrons 
are treated fully relativistically and the valence electrons 
scalar-relativistically. The LAPW basis consists of 
atomic-like functions within MT spheres at the atomic 
sites and plane waves in the interstitial region. The 
conventional basis set is furthermore expanded with 
local orbitalsS where appropriate. Inclusion of local 
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orbitals in addition to the normal FLAPW basis enforces 
mutual state orthogonality and increases variational 
freedom. 

The explicit MLWF construction was performed with 
the corresponding publicly available coda^. Several mi- 
nor additions to the exisiting code were performed in 
this work in order to account for the specifc interfacing 
requirements within LDA+DMFT. 

The NMTO construction was performed on the basis 
of scalar-relativistic LMTOS calculations in the atomic- 
sphere approximation (ASA) with combined corrections. 
Also LMTO is an all-electron method, i.e., it is fully self- 
consistent for core and valence electrons. We utilized the 
Stuttgart TB-LMTO-ASA codeS^. 

2,. Impurity-model solver 

The crucial part of the DMFT framework is the solu- 
tion of the effective quantum impurity problem. Depend- 
ing on the symmetries of the specific case at hand, and 
the demands for accuracy, several different techniques are 
available to solve this problem in practice (for reviews 
see Ref. |6l9j |'). First the on-site interaction vertex has 
to be defined. In both cases, i.e., SrVOa and BaVSs, we 
are facing a realistic three-band problem. We keep only 
density-density interactions in Hu, thus no spin- flip or 
pair-hopping terms are included. When neglecting ex- 
plicit orbital dependence of the interaction integrals, Hu 
reads then as 

U' 

Hu = U^flrn-fflnil + — ^ flmaflra' S 

rn mm' a 

m^m' 

U" 

+ — ^ flmafljn'a- ■ (26) 

mim'a 
m^m' 

Here hjn(j=dlj^„d^^, where m,a denote orbital and spin 
index. The following parametrization of U' and U" has 
been proven to be reliable^SiSl in the case of i2g-based 
systems: U'=U—2J and U"=U—3J. No explicit double- 
counting term Sdc was introduced in our specific calcula- 
tions. This is due to the fact that we used C=>V, i.e., our 
correlation subspace was chosen to be identical with the 
set of Wannier bands. In that case the double counting 
may be absorbed in the overall chemical potential. 

The solution of the quantum impurity problem corre- 
sponds to the evaluation of the impurity Green's func- 
tion Giijip for a given input of the dynamical mean-field 
(Eq. H1U|) '). which may be expressed within the path- 
integral formalism via 

~ r') ^ -(r4.(r)4(T'))5.„,, (27) 
with Zi„p= fY[V{d^;d} e"^'-^ , 



where S'imp is the effective action defined in Eq. ||HJ|. We 
utilize the auxiliary-field Quantum Monte Carlo (QMC) 
method following Hirsch-Fye^S to compute H27|l . In this 
method the path integral is evaluated by a stochastic in- 
tegration. Therefore S'imp is represented on L discrctized 
imaginary time slices of size Ar=/?/L. Since the vertex 
Hu is quartic in the fermionic degrees of freedom, a de- 
coupling using an exact discrete Hubbard-Stratonovich 
transformation is needed. For M orbitals involved, a 
number of AI{2M — 1) so-called "Ising fields" emerge 
from this decoupling for each time slice. In the end, the 
number of time slices L and the number of Monte-Carlo 
sweeps iVivic arc the sole convergence parameters of the 
problem. The QMC technique has no formal approxima- 
tions, however the numerical effort scales badly with M 
and p. 

Note again that although we so far outlined 
LDA-I-DMFT as a fully self-consistent scheme, i.e., in- 
cluding charge-density updates, the results in the follow- 
ing sections were obtained from a simpler post-processing 
approach. Thereby the self-consistent LDA Wannier 
Hamiltonian was used in H22|) and no charge-density up- 
dates were performed. 

III. RESULTS 
A. SrVOa 

1. Characterization and band- structure calculations 

A quadratic temperature behavior of the resistivity 
up to room temperature^^, albeit with a large pref- 
actor. qualifies the electronic structure of the 3(i(t2g)^ 
compound SrVOs as a Fermi liquid with intermediate 
strength of the electron-electron interactions. Still, a di- 
rect comparison of the photoemission spectral function 
with the one-particle density of states (DOS), calculated 
e.g. within DFT-LDA, yields poor agreement, indicat- 
ing a strong need for an explicit many-body treatment of 
correlations effects. DFT-LDA also yields a specific-heat 
coefficient (slope of C/T at low-T) which is too small by 
approximately a factor of two, i.e., the electronic effective 
mass is enhanced due to correlation effects. 

A perfectly cubic perovskite structure and the ab- 
sence of magnetic ordering down to low temperatures 
makes SrVOs an ideal test material for first-principles 
many-body techniques. It has thus been the subject of 
many experimental and theoretical investigations (using 
IDA [ pjiPT )^^-^^-'^^-^''-'^^-^^-''^-''^-''^-''^-''^-''^-''^-'''' In the 
present work, we use this material as a benchmark for 
our Wannier implementation of LDA-I-DMFT. with re- 
sults very similar to previous theoretical studies. 

We start the investigation of SrVOs with a brief DFT- 
LDA study. The crystal structure of the transition-metal 
oxide SrVOa is rather simple, exhibiting full cubic sym- 
metry (space group Pm3m) with a measuredS^ lattice 
constant of 7.2605 a.u.. The V ion is placed in the center 
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FIG. 2: (Color online) (a) LDA data for SrVOs calculated 
with the MBPP code, (a) Band structure, (b) DOS. For the 
local V(3ci) /0(2p)-D0S the cutoff radius was half the nearest- 
neighbor distance, respectively. 



of an ideal octahedron formed by the surrounding O ions 
see Fig. 0] The O ions are at the face centers of a cube 
having V at its center and Sr at its corners. 

Figurc|21shows the band structure and the DOS within 
LDA. The data reveals that there is an isolated group of 
partially occupied bands at the Fermi level, with a total 
bandwidth of 2.5 eV. For an ion at site R the local orbital 
density matrix n|jjt,/~I]ki/ /ki'(V'kiy|-f^^m IV'ki.) is a 

measure for the occupation probabilities within the set 

of, say cubic, harmonics {Ki\^}. In the case of SrVO^ 
this matrix is diagonal, and it is seen in Fig. ^jp that 
from such a projection the bands at £p may be described 
as stemming dominantly from V(t2g) orbitals. Since 
the three t2g orbitals, i.e, {d^y, dxz, dyz} are degenerate, 
they have equal contribution to the bands. Due to the 
full cubic symmetry the distinct t2g orbitals are nearly 
exclusively restricted to perpendicular planes which ex- 
plains the prominent 2D-like logarithmic-peak shape of 
the DOS. The V(eg) states have major weight above the 
Fermi level, whereas the 0(2p) states dominantly form 
a block of bands below £p . The energy gap between 
the 0{2p) and t2g block amounts to 1.1 eV. In spite of 
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FIG. 3: (Color online) t2g bands for SrVOa using dif- 
ferent schemes to compute the t2g Wannier functions 
(and the underlying LDA band structure or poten- 
tial), dark: MLWF(MBPP), dashed-red (dashed-gray) ; 
MLWF(FLAPW) and green (light gray): NMTO(LMTO- 
ASA). The t2g bandwidth is marginally larger in FLAPW, 
leading to small differences. 

the "block" characterization, there is still significant hy- 
bridization between the most relevant orbitals, i.e., V(3c?) 
and 0(2p), over a broad energy range. 



2. Wannier functions 

The low-energy physics of SrVOa is mainly determined 
by the isolated set of three t2g-like bands around the 
Fermi level. This suggests the construction of an effective 
three-band Wannier Hamiltonian as the relevant minimal 
low-energy model. In the following, we construct Wan- 
nier functions associated with this group of bands, and 
also pick these three Wannier functions as generating the 
correlated subset C, so that W=C in the notations of the 
previous section. This choice of course implies that the 
resulting Wannier functions, though centered on a vana- 
dium site, have also significant weight on neighbouring 
oxygen sites. More localized functions can indeed be ob- 
tained by keeping more bands in the Wannier construc- 
tions (i.e. by enlarging the energy window) and thus 
keeping W larger than C, as described at the end of this 
subsection. However, we choose here to explore this min- 
imal construction as a basis for a DMFT treatment and 
show that it actually gives a reasonable description of 
this material. 

Figure |3| exhibits the Wannier bands obtained within 
our three utilized schemes: maximally-localized WFs 
from the MBPP and FLAPW codes (abbreviated 
in the following respectively by MLWF(MBPP) and 
MLWF (FLAPW)) and the NMTO scheme used as a post- 
processing tool on top of the LMTO-ASA code (denoted 
as NMTO(LMTO-ASA)). For the MLWF construction 
a starting guess for the WFs was provided by utilizing 
atomic-like functions with t2g symmetry centered on the 
V site. Some details on the construction of the NMTO- 
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FIG. 4: (Color online) faj-like MLWF w^y for SrVOa de- 
rived from the MBPP code. First row; SrVOa structure with 
Sr (large blue/dark), V (red/gray) and O (small yellow/light 
gray) and perspective view on Wxy Second row: Wxy viewed 
along the c axis and along a axis. The contour value w^J was 
chosen as 0.05 (a.u.)"''''^. 

WFs are provided in Appendix [CI Both MLWF and 
NMTO schemes yield bands identical to the LDA bands. 
The small discrepancies seen in Fig. 13 are due to differ- 
ences in the self-consistent LDA potentials. This overall 
agreement between the different methods reflects the co- 
herent LDA description for this material. 

Although all three sets of WFs span the same Hilbcrt 
space, and the bands are therefore the same, the MLWFs 
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FIG. 5: (Color online) Contour-lines plot of the t2g-\ike 
MLWF Wxy for SrVOs. Distinct contour values (in (a.u.)"^/^) 
are given in the plot. 
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FIG. 6: (Color online) (a) t29-like w^y WF. (a) From V to V 
along [110], and (b) from O to O along [010]. 



and the WFs obtained by symmetrically orthonormaliz- 
ing the NMTO set are not necessarily identical. In order 
to compare the Wannier orbitals, we generated the set 
{w„(r)} within a (3x3x3) supercell on a (120x120x120) 
real-space mesh. As an example, Fig.01shows the t2g-\ike 

Wxy{v) Wannier orbital for a chosen constant value w^xy 
as obtained from the MLWF(MBPP) construction. By 
symmetry, all three Wannier orbitals come out to be cen- 
tered on the V site. A general contour plot for Wxy{v) is 
given in Fig. |S1 The Wannier orbitals show clear t2g sym- 
metry, but in addition have substantial oxygen character, 
7r-0(2p) in particular. The important hybridization be- 
tween the V(t2g) and 0{2p) atomic- like orbitals seen in 
Fig. [2| is explicitly transfered in the Wannier orbital. By 
comparing the three different sets of Wannier orbitals we 
find remarkably close agreement. Thus the MLWF and 
NMTO constructions provide nearly identical vanadium 
t2g Wannier orbitals in the case of cubic SrVOs. A de- 
tailed comparison is shown in Fig. where the WFs are 
plotted along specific directions. From these graphs it 
may be seen that the MLWF (MBPP) shghtly disagrees 
with the WFs from the two other schemes close to the 
nuclei. This discrepancy is due to the pseudization of the 
crystal wave functions close to the nucleus. Although the 
3(i wave function is nodeless, the pseudo wave function 
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TABLE I: Spread (r^} of the V-centered tag WFs for SrVOs. 
We employed a (8x8x8) fc-point mesh for the construction 
of the MLWFs. The spreads were calculated in two ways. 
First from the fc-space integration performed in in the MLWF 
code, and second via an r-space integration within a (3x3x3) 
supercell. Additionally shown is the respective normalization 
of the WFs within this supercell. Finally, for completeness 
the Q value minimized within the MLWF construction is also 
given, even if due to the complete degeneracy it reduces to 
three times the spread (r^). 



scheme 


(r^) (a.u.2) 


norm 


n (a.u.^) 


k space r space 


r space 


k space 


MLWF(MBPP) 

MLWF(FLAPW) 

NMTO(LMTO-ASA) 


6.86 6.64 
6.96 6.75 
6.82 


0.998 
0.997 
0.995 


20.57 
20.93 


Ref. ^ 


8.46 







is modified in order to provide an optimized normcon- 
serving pseudopotential. However, this difference in the 
WFs has no observable effect on the description of the 
bonding properties as outlined in general pseudopoten- 
tial theorji^Sia (see also Tab.|llll. Only marginal differ- 
ences between the different WFs can be observed away 
from the nuclei. Generally, the fast decay of the WFs is 
documented in Fig. In this respect, Table exhibits 
the values for the spread (r^) of the WFs from the differ- 
ent schemes. The MBPP and FLAPW implementations 
of the MLWFs have spreads which differ by 2%. Since 
for the MLWFs the spread has been minimized, that of 
the NMTO-WFs should be larger, and it indeed is, but 
merely by a few per cent. So in this case the NMTO-WFs 
may be seen as maximally localized, also in the sense of 
Ref. (2^. A substantially larger value for the spread is 
however obtained from the orthonormal LMTOs, as seen 
from Ref. ,22] . 

To finally conclude this part of the comparison, we de- 
duced the relevant near-neighbor hopping integrals from 
the real-space Hamiltonian in the respective Wannier ba- 
sis, given in Table. |n] The dominance of the nearest- 
neighbor hopping in connection with the fast decay of 
the remaining hoppings clearly demonstrates the strong 
short-range bonding in SrVOs. The close agreement 
of the hoppings between the three different Wannier 
schemes again underlines their coherent description of 
this material. It can be concluded that although con- 
ceptually rather different. MLWF and NMTO provide a 
nearly identical minimal Wannier description for SrVOa. 
The small numerical differences seem to stem mainly 
from the differences in the electronic-structure descrip- 
tion within the distinct band-structure methods. 

At the end of this subsection we want to draw atten- 
tion to the fact that the performed minimal Wannier con- 
struction solely for the t2g bands is of course not the only 
one possible, as already mentioned above. Depicted in 
Fig. [3 are the WFs obtained by downfolding the LDA 
electronic structure of SrVOa to V{3d) and 0{2p) states. 
Hence this corresponds to describing SrVOs via an 14- 



TABLE IL Symmetry-inequivalent intersite hopping integrals 
-ffyz.j/z for SrVOa. Energies are in meV. 



xyz 


001 


100 


Oil 


101 


111 


002 


200 


MLWF(MBPP) 


-260.5 


-28.2 


-83.1 


6.5 


-6.0 


8.4 


0.1 


MLWF(FLAPW) 


-266.8 


-29.2 


-87.6 


6.4 


-6.1 


8.3 


0.1 


NMTO(LMTO-ASA) 


-264.6 


-27.2 


-84.4 


7.3 


-7.6 


12.9 


3.5 




FIG. 7: (Color online) Distinct WFs for SrVOa obtained 
from the MLWF construction using the MBPP code. First 
row: 0{px), 0{py) and 0(pz) for a chosen oxygen site. Sec- 
ond row: V{t2g,xy) as well as V{eg,3z^-r^) and Y{eg,x^-y^). 
The contour value for each of the MLWFs was chosen as 0.05 
(a.u.)-3/^ 



band model, i.e., three p orbitals for three O ions and 
five d orbitals for the single V ion in the unit cell. Due to 
minor degeneracies with higher lying bands (see Fig. ^ 
the disentangling procedure for the MLWF construction 
has to be used, but no relevant impact is detected in this 
case. Now there are distinct WFs for 0(2p) and V(3d) 
with significantly smaller spreads. Individually the lat- 
ter are in a.u.^: 2.61 for V(t2g) and 2.32 for V(eg). and 
2.68 for cr-0(2p) and 3.39 for 7r-0(2p), resulting in a total 
spread of fi=40.75 a.u.^. 



3. LDA +DMFT calculations 

Thanks to the simplicity of the perfectly cubic per- 
ovskite structure and the resulting degeneracy of the 
three t2g orbitals, SrVOa is a simple testing ground for 
first-principles dynamical mean-field techniques. In fact 
SrVOs is quite a unique case in which the calculation of 
the local Green's function H22|l . which usually involves 
a k summation, can be reduced to the simpler calcula- 
tion of a Hilbert transform of the LDA DOS. Indeed, 
because of the perfect cubic symmetry, all local quan- 
tities in the t2g subspace are proportional to the unit 
matrix: Gj5^^,(iw„) GiociiuJn)Smm' , AE„™'(ia;„) = 
AT,{iuJn)Sm7n' (as wcll as the LDA DOS D,nm'{£) ~ 
D{e)Smm' projected onto the orbitals Xm), so that H22() 
reduces to: Gioc{i^n)—J +i!^e-AY. - Note however 
that this does not hold in general for other materials, 
as soon as the local quantities are no longer proportional 
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to the unit matrix. Although many actual LDA+DMFT 
calculations in literature use this representation as an ap- 
proximation to the correct form given by Eq. In the 
calculations documented in this work we always used the 
more generic Hamiltonian representation and k summa- 
tions. 

Taking into account the strong correlations within the 
t2g manifold results in substantial changes of the local 
spectral function compared to the LDA DOS, namely 
a narrowing of the QP bands close to the Fermi level 
while the remaining spectral weight is shifted to Hub- 
bard bands at higher energies. This general physical pic- 
ture of the correlated metal can be understood already 
in the framework of the multi-orbital Hubbard model 
as the coexistence of QP bands with atomic-like exci- 
tations at higher energy. It directly carries through to 
the realistic case of SrVOa as studied in several previous 
— crln^^-^''''^^-^^. 

Moreover, an important feature of LDA+DMFT that 
emerges in the present case of a completely orbitally de- 
generate self-energy has been put to test against experi- 
ments. Indeed, in this special case Fermi-liquid behavior 
in conjunction with a k-indcpcndcnt self-energy leads to 
the value of the local spectral function p(e) at the Fermi 
level being equal to its non-interacting counterpart just 
in the same way as in the one-band Hubbard modetf^. 

In this work we performed LDA-I-DMFT calcula- 
tions for SrVOs by using the self-consistent LDA Wan- 
nicr Hamiltonian Hks derived from the different band- 
structure codes, i.e., MBPP, LMTO-ASA and FLAPW, 
described above. As expected from the good agreement 
of the band structure and the Hamiltonians the result- 
ing Green's functions are identical within the statistical 
errors bars (see the inset of Fig. IHJ. Fig. IHlalso displays 
the local spectral functions based on the MLWF(MBPP) 
scheme and calculated for different values of U. The "pin- 
ning" of p(0), independently of the value of the inter- 
actions is clearly visible, despite the finite temperature 
of the calculations. This indicates that the calculations 
have indeed been performed at a temperature smaller 
than the QP (Fermi-liquid) coherence scale of this mate- 
rial. 

Figurcinidisplays the local spectral function convoluted 
with an assumed experimental resolution of 0.15 eV and 
multiplied by the Fermi function. This quantity repre- 
sents thus a direct comparison to angle-integrated pho- 
toemission spectra (albeit neglecting matrix elements, 
which can in certain circumstances appreciably depend 
e.g. on the polarization of the photons, see Ref. |8(j|'). 

The general agreement with recent experimental 
data22iSi2Si2L2£ is reasonable. Photoemission experi- 
ments locate the lower and upper Hubbard bands at ener- 
gies about -2 eV to -1.5 e V'^''i^° and 2.5 eV— respectively. 
In our calculations the lower Hubbard band extends be- 
tween -2 eV to -1.5 eV, while the maximum of the up- 
per Hubbard band is located at about 2.5 eV, for values 
of the Coulomb interaction U of about 4 eV. However, 
we also confirm the findings of Ref. [sOl who point out 
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FIG. 8: (Color online) Spectral function for SrV03 for 
different values of U resulting from LDA+DMFT in the 
MBPP implementation. Inset: comparison of Green's func- 
tions from MBPP, FLAPW and NMTO implementations of 
LDA+DMFT for (7=4 eV. In all calculations J=0.65 eV was 
used. 



that LDA+DMFT calculations generally locate the lower 
Hubbard band at slightly higher (in absolute value) bind- 
ing energies than -1.5 eV, the energy where their data 
exhibits its maximum. 

Concerning the choice of the Coulomb interaction U 
different points of view can be adopted. First, one can 
of course choose to try to calculate U itself from first 
principles by e.g. constrained LDA224Si21S2& or RPA- 
based techniquesSLSS. Another option is to use it as an 
adjustable parameter and to determine it thus indirectly 
from experiments. While in the present case the order 
of magnitude of the interaction (t/^3.5-5.5 eV)S2i2LSS is 
indeed known from first-principles approaches, the exact 
values determined from different methods still present a 
too large spread to be satisfactory for precise quantita- 
tive predictions. We therefore adopt the second point of 
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FIG. 9: SrVOa spectral function convoluted with a Gaussian 
experimental resolution (assumed to be 0.15 eV) and with the 
Fermi function at /3=30 eV^^. 
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view here, noting that U values of around 4 eV repro- 
duce well the experimentally observed2Si2Si2£ mass en- 
hancement of ^1.8 to 2. The agreement concerning the 
position of the Hubbard bands seems to be fair, given the 
theoretical uncertainty linked to the analytical continua- 
tion procedure by maximum-entropy techniques and the 
spread in available experimental data. Still, it is conceiv- 
able that the determination of the precise position of the 
Hubbard bands could require more sophisticated meth- 
ods than LDA-I-DMFT done with a static U parameter, 
and that in fact we arc facing the consequences of sub- 
tle screening effects which, within an effective three-band 
model, could only be described by a frequency-dependent 
interactional. 



B. Cmc2i-BaVS3 

1. Structure and physical properties 

The transition- metal sulfide BaVSa is also a 3rf(i2g)^ 
system, but its physical properties are far more com- 
pjgjj9ij92 \\^QT^ those for the cubic perovskite SrVOa con- 
sidered above. In a recent work22i24, three of us have sug- 
gested that a correlation-induced redistribution of orbital 
populations is the key mechanism making the transition 
into a charge-density wave (CDW) insulating phase pos- 
sible. Here, we use our Wannier formalism to make a 
much more refined study of this phenomenon and to cal- 




FIG. 10: (Color online) BaVSs in the Cmc2i structure. The 
V ions are shown as smaller (red/gray) spheres, the Ba ions 
as larger (blue/dark) spheres. 



culate how correlations modify the Fermi-surface sheets 
of the metal. Also, this material is a challenging test- 
ing ground for the Wannier construction because of the 
strong hybridization between the transition-metal and 
ligand bands. 

We first give a very brief summary of some of the phys- 
ical properties of BaVSs of relevance to the present pa- 
per. At room temperature BaVSa exists in a hexagonal 
crystal structure (space group P63 / mmc) , with two for- 
mula units of BaVSa in the primitive cell. There are 
straight chains of face sharing VSg octahedra along the 
c axis, and Ba ions in between. A continuous structural 
phase transition at Ts'^240 K reduces the crystal sym- 
metry to orthorhombic, thereby stabilizing the Cmc2i 
structure, again with two formula units in the primitive 
cell. Now the VS3 chains are zigzag-distorted in the be 
plane. In this phase. BaVSs is a quite bad metal, with 
unusual properties such as a Curie- Weiss susceptibility 
from which the presence of local moments can be inferred. 
At ~70 K a second continuous phase transition takes 
place2SiS&, namely a metal-insulator transition (MIT) be- 
low which BaVSs becomes a paramagnetic insulator. A 
doubling of the primitive unit cell22i22iSS is accompanying 
the MIT. Together with large one-dimensional structural 
fluctuations along the chainsS^ and additional precursive 
behavior for the Hall constanliSi just above Tmit, the 
transition scenario is reminescent of a Peierls transition 
into a charge-density wave (CDW) state. Finally a third 
second order transition appears to occur at ~30 K. This 
so-called "X" transition is of magnetic kind and shall an- 
nounce the onset of incommensurate antiferromagnetic 
orderiiffi in the insulator. 

Here we want to focus on the orthorhombic (Cmc2i) 
structure (see Fig. IIUII at r=100 K, i.e., just above the 
MIT. Ten ions are incorporated in the primitive cell. 
Whereas the two Ba and two V ions occupy (4a) sites, 
there are two types of sulfur ions. Two S(l) ions are po- 
sitioned at (4a) apical sites on the b axis, while four S(2) 
ions occupy (8b) sites. The lattice parameters are>i2i 
a=12.7693 a.u., 6=21.7065 a.u. and c=10.5813 a.u.. 



2. Band structure 

Figure ^2 depicts the LDA band structure and DOS 
of Cmc2i-BaVS3. To allow for orbital resolution, the 
local DOS was again projected onto symmetry-adapted 
cubic harmonics by diagonalizing n^^,. It is seen that 

^ ^ mm 

the bands at have dominant t2g character, however 
they still carry sizeable S(3p) weight. Furthermore, the 
t2g-like bands are now not isolated but strongly entangled 
with S(3p)-like bands. Due to the reduction of symmetry 
from hexagonal to orthorhombic, the t2g manifold splits 
into Aig and Eg . The two distinct Eg states will be 
denoted in the following Egi and Eg2. Being directed 
along the c axis, the Aig orbital points towards neigbor- 
ing V ions within a chain and the corresponding band 
(see Fig. IT^ ) shows a folded structure because of the 
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FIG. 11: (Color online) LDA data for BaVSs after the 
MBPP method, (a) band structure and(b) DOS. For the 
local V(3d)/S(3p)-DOS the cutoff radius was half the min- 
imum nearest-neighbor distance, respectively. The inset in 
(b) shows the symmetry-adapted local V(3d)-DOS close to 
the Fermi level. 



existence of two symmetry-equivalent V ions in the unit 
cell. The folded Aig band has a bandwith of 2.7 eV, while 
Eg states form very narrow (0.66 cV) bands right at the 
Fermi level. 

From the LDA DOS it seems that a projection onto 
t2g-like orbitals close to the Fermi level by diagonalizing 

"-mm' least to a first approximation, meaningful^Si. 

However there is a substantial S(3p) contribution close 
to s-p and generally large charge contributions in the in- 
terstitial. Hence establishing a very accurate correspon- 
dance between relevant bands and orbitals is not pos- 
sible in such a way. In contrast, the Wannier schemes 
discussed above are quite suitable for dealing with this 
situation. To be specific, we applied the MLWF scheme 
in the new MBPP implementation to this problem. 

Besides providing a test for the MLWF scheme, the 
present study will allow us to make considerably more 
precise the findings of Ref. [o^ regarding the crucial role 
of correlation-induced changes in the orbital populations, 
and most notably to clarify how these changes can mod- 



ify the Fermi surface (FS) of this material in such a way 
that favorable conditions for a CDW transition indeed 
hold. Key to the physics of BaVSs is the simultane- 
ous presence of two quite distinct low-energy states, the 
rather delocalized Aig and quite localized Eg , among 
which the electronic density with one electron per vana- 
dium has to divide itself. Depending on temperature, 
the associated orbital populations correspond to the best 
compromise between gain of kinetic energy and cost of 
potential energy. As it appears, this compromise seems 
to be realized by a CDW state below the MIT. However, 
as revealed in several electronic structure studiesSSiSiiS^, 
a DFT-LDA description of BaVSs does not explain the 
occurence of a CDW instability. Though the mainly Aig- 
like band appears to be a promising candidate, a nest- 
ing scenario in agreement with the critical wave vector 
qc=0.5c* from experiment^ is not realizable. In Fig. 112b 
we elaborated a so-called "fatband"ii2^ resolution of the 
LDA band structure close to the Fermi level, which is 
helpful to reveal the respective band character to a good 
approximation. Thereby the Bloch function associated 
with a given /j-point and eigenvalue is projected onto or- 
thonormal symmetry-adapted local orbitals (determined 
as usual by diagonalizing the local orbital density matrix 
^mm')- T^^c resulting magnitude of the overlap is de- 
picted as a broadening of the corresponding band. Here 
it is seen that Aig-like band cuts the Fermi level close 
to the boundary of the BZ along F-Z, i.e., along the c* 
axis in reciprocal space. Since the Z point is located at 
kz=0.5c*, in numbers this amounts to 2kp~0.94: for the 
Aig-like band within LDA, nearly twice the experimen- 
tal value determined for the nesting vector. Furthermore, 
also other parts of the LDA FS are out of reach for qc, as 
the Aig sheet is too extended and additionally strongly 
warped (see also Fig. Iiyb). In other words, LDA appar- 
ently overestimates the population of the more itinerant 
Aig state. 



Moreover the role of the electrons with strong Eg char- 
acter at the MIT is not obvious. When approaching 
Tmit these nearly localized electrons should surely con- 
tribute to the Curie- Weiss form of the magnetic sus- 
ceptibility2i2&. In fact the "bad-metal" regimeS^ above 
the MIT, including significant changes in the Hall coeffi- 
cientSi, might largely originate from scattering processes 
involving the Eg electrons. But even if the Aig bands be- 
come gapped at the MIT, from an effective single-particle 
LDA viewpoint the remaining Eg bands may still en- 
sure the metallicity of the system. We therefore believe 
for several reasons that correlation effects beyond LDA 
are important^^ for an understanding of the physics of 
BaVSa. We will further outline relevant mechanisms, 
now based on a more elaborate Wannier scheme, in sec- 
tion |mB4l 
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FIG. 12: (Color online) (a) t2g fatband resolved band struc- 
ture of BaVSa. The color code is as follows: Aig (blue/dark), 
Egi (red/gray) and £32 (green/light gray), (b) downfolded 
t2g Wannier bands (light blue/light gray) for BaVSa obtained 
from the MLWF construction using the MBPP code. 



3. Wannier functions 

The central difRculty in constructing t2g-\ike WFs for 
BaVSa is the strong hybridization between V{3d) and 
S(3p), leading to a substantial entanglement between 
the two band manifolds. In detail, whereas the two Eg 
states form four very narrow bands, mainly confined to 
the Fermi level, the folded Aig band extends into the 
dominantly S(3p)/V(eg) band manifolds lower/higher in 
energy. This entanglement is documented in Fig. 112b 
by significant "jumps" of the corresponding Aig fatband 
between different bands. One may of course downfold 
the BaVSa band structure including not only V{t2g) but 
also S(3p) and V(eg) orbitals. However, in this work we 
wanted to investigate the properties and reliability of the 
minimal, i.e., i2g-only, model. In the following we discuss 
the results obtained via the MLWF construction. Cor- 
responding studies were also performed with an NMTO 
basis set leading to the same physical picture. But a de- 
tailed comparison would at this point shift the attention 
from the investigated physical mechanisms. 

In order to downfold onto {Aig, Egi, £"^2} we employed 
the disentangling proceduro^ of the MLWF construc- 
tion. The WFs were initialized via cubic harmonics 
adapted to an ideal local hexagonal symmetry. To the 
aim of correct disentangling of the six Wannier target 
bands we provided twenty bands in an outer energy 
window around the Fermi level for the construction of 
Mj^J^;*'^ ■ III order to reproduce the LDA FS and the band 



dispersions close to the Fermi level correctly, we addi- 
tionally forced the Wannier bands in an inner energy 
window near e-p to coincide with the true LDA bands^. 
The initial WFs correspond to an optimized f2/=101.58 
a.u.^ and a starting value f2=3.62 a.u.^, hence a total 

of 105.19 a.u.^ for the chosen energy windows. Af- 
ter '^50000 iteration steps finally converged to 103.30 
a.u.^. During the minimization process, adaptation of 
the WFs to the true orthorhombic symmetry was clearly 
observed by the occurence of distinct steps in il. The 
resulting t2g Wannier bands are shown in Fig. 112b in 
comparison with the original LDA band structure. It 
is seen that the Wannier bands at the Fermi level are 
truly pinned to the original LDA bands. Furthermore the 
interpolated lowest/highest Wannier band follows nicely 
the former Aig fatbands. The same t2g dispersion is also 
obtained within an NMTO contruction. 

Table IITTI reveals that the Wannier spreads are signifi- 
cantly larger than for the t2g WFs in SrVOa. The spread 
for the Aig WF is slightly smaller than those for the 
Eg orbitals. Moreover, it is seen from Tab. IIIII that the 
Wannier centers are shifted from the V sites. This may 
be explained by the low symmetry of the Cmc2i-BaVS3 
structure, since already from the missing inversion sym- 
metry in the site symmetry for the V ion there is no need 
for the WFs to be centered on the V sites^S. Additionally, 
the explicit inclusion of the sulfur contribution to the t2g 
WFs leads to directed Wannier orbitals. The latter may 
be inspected in Fig. ^1 where we plotted the three Wan- 
nier orbitals associated with the first V ion in the unit 
cell (the other three are directly related by symmetry due 
to the equivalence of the V ions). Indeed the Aig-like 
orbital is directed along the c axis, whereby the zigzag 
distortion of the chains causes some tilting. Note that 
the orbitals have some weight on neighboring V sites. As 
especially observable for the E'gi-like orbital, this weight 



TABLE III: Wannier centers Ru, and spread (r^) of t2g- 
like MLWFs for BaVSs constructed from a (6x6x6) fc-point 
mesh. The positions of the symmetrically equivalent V 
sites in cartesian coordinates read Rv(i) = (0. 00, 0.46,-0. 01) 
a.u. and Rv(2) =(0.00,-0.46,5.28) a.u.. The V(2) site is 
symmetry-related to the V(l) site by the symmetry opera- 
tion C;^^'Rv(i)+0.5. 



WF 


Ru, (a.u.) 


Hw — Rv (a.u.) 


(r'} (a.u.2) 


Aig, V(l) 
Egl, V(l) 
Eg2, V(l) 


0.00, 0.75, -0.20 
0.00, 0.64, 0.38 
0.00, 1.02, -0.32 


0.00, 0.30, -0.19 
0.00, 0.18, 0.39 
0.00, 0.56, -0.31 


16.57 
17.55 
17.53 


TABLE IV: Wannier centers Ru, and spread {r^) as in 
Tab, lini now in the crystal-field basis. 


WF 


Ru, (a.u.) 


Hw — Rv (a.u.) 


(r') (a.u.^) 


Aig, V(l) 
Egl, V(l) 
Eg2, V(l) 


0.00, 0.75, -0.17 
0.00, 0.65, 0.34 
0.00, 1.02, -0.32 


0.00, 0.30, -0.16 
0.00, 0.19, 0.35 
0.00, 0.56, -0.31 


16.60 
17.55 
17.53 
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FIG. 13: (Color online) t2g Wannier orbitals for BaVSs form- 
ing the maximally-localized basis. The columns from left to 
right show the Aig, Egi and Eg2 orbitals, while the second 
row displays the orbitals viewed along the c axis. There the a 
axis is vertically oriented, while the b axis horizontally. The 
contour value for each of the orbitals was chosen as 0.045 
(a.u.)-3/^ 







FIG. 14: (Color online) t2g Wannier orbitals for BaVSa form- 
ing the crystal-field basis. The order of the columns and row 
same as in Fie ll3l Note that the Eg2 orbital remains invari- 
ant under this transformation. The contour value for each of 
the orbitals was equally chosen as 0.045 (a.u.)"^'^^. 



nects to the apical S(l). For Aig the sulfur distribution 
is in-between, yet favoring the nearest-neighbor S(l). 

In order to investigate the properties of these WFs in 
more detail. Tab. displays the hopping integrals for 
relevant paths on the lattice. One immediately realizes 
that while the Eg2 orbital remains isolated, there is a 
sizable hybridization between the Aig and Egi Wannier 
orbitals. The latter hybridization clearly couples Aig and 
Egi in view of the zigzag distortion along the b axis (see 
Fig. I13|l This coupling should play an important role for 
the understanding of the whole BaVSa phase diagram, 
also crucial for the hexagonal-to-orthorhombic structural 
transitioniSi. Note also the dominant Egi-Egi hopping 
along the a axis, i.e. [110], which is in accordance with 
orbital extension of the corresponding WFs (see Fig. ll3|l . 

In Figure [T^ we plot the disentangled (effective) t2g 
bands with a fatband resolution to reveal the respective 
contribution of the obtained MLWFs (for simplicity we 
do not plot the Eg2 fatband). One can see that the Aig 
band carries substantial Egi weight (and vice versa) due 
to their sizeable hybridization in the maximally-localized 
basis. This is is somehow counter-intuitive to the original 
low-energy picture of a broader Aig band and narrower 
Egi band. Though some minor Aig-Egi hybridization 
in the F-C-Y plane is in line with the original LDA fat- 
bands, it seems that this hybridization is slightly over- 
represented in the maximally-localized basis. Surely, the 
latter basis is not a priori physically designated, and it 
might be that the the straightforward MLWF construc- 
tion does not provide the most suitable physical Wannier 
basis for the description of BaVSs. 

With the aim of reducing this strong Aig-Egi 



TABLE V: Hopping integrals between the t2g Wannier or- 
bitals of BaVSs in the maximally-localized basis. The term 
'00 1' shall denote the hopping to the nearest-neighbor V site 
within the unit cell. One of the nearest-neighbor V ions in the 
ab plane is located at '100', while '110' and 'IlO' are closest 
V ions along a and b, respectively. Energies in moV. 





Alg-Alg 


Egl-Egl 


Eg2-Eg2 


Alg-Egl 


Alg-Eg2 


Egl-Eg2 


000 


AUA 


218.0 


235.6 


40.8 


0.0 


0.0 


00 i 


-441.5 


-24.7 


-12.4 


-242.6 


0.0 


0.0 


001 


-66.0 


-5.4 


2.7 


-8.9 


0.0 


0.0 


100 


-30.4 


8.5 


-26.1 


-15.8 


16.4 


-10.9 


no 


-17.4 


-84.4 


29.2 


-16.5 


-0.7 


11.6 


110 


1.6 


1.7 


-5.5 


-1.8 


0.0 


0.0 



has local eg symmetry. Thus the MLWF construction re- 
produces here the intuitively formulated symmetry con- 
straints imposed in the NMTO construction2i*^2i2^. This 
symmetry relation between WFs on neigboring sites has 
also been noted for the V2O3 compouncii24, half- metallic 
ferromagnet9i2^ and, most dramatically, for Na2:Co02 
compoundsiS^. The Egi orbital is mainly oriented in 
the plane defined by the S(2) ions, while Eg2 remains in 
the corresponding perpendicular plane. Accordingly, Egi 
has stronger weight on the S(2) ions, whereas Eg2 con- 



TABLE VI: Hopping integrals as in Tab. IVl but now in the 
crystal-field basis. 





Alg-Alg 


Egl-Egl 


Eg2-Eg2 


Alg-Egl 


Alg-Eg2 


Egl-Eg2 


000 


422.6 


209.8 


235.6 


0.0 


0.0 


0.0 


00 i 


-510.5 


44.3 


-12.4 


-145.7 


0.0 


0.0 


001 


-85.5 


14.2 


2.7 


7.1 


0.0 


0.0 


100 


-35.4 


13.6 


-26.1 


-7.0 


-14.0 


13.9 


110 


-26.3 


-75.5 


29.2 


-28.1 


-1.6 


-11.5 


110 


1.2 


2.1 


-5.5 


-1.7 


0.0 


0.0 
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FIG. 15; (Color online) t2g Wannier fatbands for BaVSs, (a) 
with respect to the maximally-localized basis, and (b) with 
respect to the crystal- field basis (see text). Color coding: 
Aig-like WF (blue/dark), £gi-like WF (red/gray). 



hybridization in the maximaUy-locahzed basis, 
we diagonahzed the on-site Wannier Hamiltonian 



H 



(T=0) 
KS 



=^jj.HKs(k) and transformed HKs(k) into 
the so-called crystal-Geld basis (e.g. also utilized in 
Ref. 01). The (again unitary) transformation is 
explicitly written as: 



(cf) 
Wp 



0.981 0.196 0.000 
-0.196 0.981 0.000 
0.000 0.000 1.000 




where the superscript 'cf marks the crystal-field basis 
and 'ml' the maximally-localized basis. This procedure 
obviously decouples Aig and Egi on average and pro- 
vides a true adaption to the local symmetry at the V 
site. Within this new basis the Aig hopping along the 
c axis is strenghtcnd at the cost of a reduced Aig-Egi 
hybridization (see Tab. IVI|I . In addition, the sign of the 
near-neighbor Egi-Egi hopping is changed from negative 
to positive. The Wannier fatbands promote now more 
the elucidated picture of the Egi bands being confined to 
Ep (see Fig. 115b ). Only minor changes may be observed 
however in the same contour plot for the transformed 
Wannier orbitals, as seen in Fig. ^1 The Aig Wannier 
orbital is now more reminescent of a d^^2_,^2 orbital per- 
pendicular to c and Egi is slightly tilted out of the plane 
defined by the S(2) ions. Correspondingly, the numbers 
for the Wannier centers and spreads have changed only 
marginally, as seen in Tab. IIVI 

By summing the spreads one easily checks that £7 is 




FIG. 16: (Color online) t2g Wannier DOS for BaVSs. (a) 
In the maximally-localized basis, and (b) in the crystal-field 
basis. 



now marginally larger than in the original maximally- 
localized basis, which is of course consistent with the fact 
that the original set was constructed by minimizing the 
spread. 

Although the real-space quantities truly do not dif- 
fer much, the electronic structure representation is very 
sensitive to rather minor changes in the basis. This is 
not only seen for the hoppings (cf. Tabs lvTVl|l but may 
also be observed when comparing the different orbitally- 
resolved DOS originating from the two Wannier basis sets 
(see Fig. I16|l . There, in general the overall lowdimen- 
sional character of the Aig band is emphasized within 
the minimal Wannier set. Note the reduced DOS magni- 
tude close to the Fermi energy for the latter band in the 
crystal-field basis. On the other hand the Aig DOS is re- 
inforced below the Fermi energy at the cost of a reduced 
Egi DOS. This effect precisely reflects the low-energy 
confinement of Egi in the crystal-field basis. 



4. LDA+DMFT calculations 

It was pointed out in Ref. |93 that correlation effects 
are important for an understanding of the MIT in BaVSa. 
From LDA+DMFT calculations using as an input the 
symmetry-adapted local DOS for {Aig^ Egi, Eg2}, a sub- 
stantial charge transfer from Aig to Eg was revealed. The 
latter was associated with a reduced Fermi wave vector 
kp for the Aig band, allowing for the possibility of the 
observed CDW instability. In this section we will check, 
and also considerably refine these earlier results. Indeed, 
because the Wannier construction allows us to study this 
intricate material within a Hamiltonian formalism which 
allows us to address k-resolved issues, we are now in a po- 
sition to study in detail the correlation-induced changes 
of the Fermi surface of BaVSa. Since we believe that 
the crystal-field basis constructed in the last section is 
more closely adapted to the physics of BaVSs than the 
direct MLWF basis, we will use this basis in the follow- 
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ing. We made cross checks using the MLWF basis, and 
wiU comment on relevant differences. 

For the new LDA+DMFT calculations we chose C/=3.5 
eV, J=0.7 eV. Although no concrete knowledge about 
the value of U for BaVSa is known, we believe that the 
latter values are in agreement with realistic values for 
this compounds^. The single-site impurity problem was 
solved for /3=30 eV'^ (~390 K). A number of 128 time 
slices and up to 10^ Monte-Carlo sweeps were used in the 
actual calculations. 

In Fig. El the orbital-resolved k-integrated spectral 
function is shown in comparision to the corresponding 
LDA DOS in the Wannier basis. The band-narrowing ef- 
fect is clearly visible. Additionally, the tendency towards 
the opening of a pseudo-like gap in the Aig band may be 
oberved. This latter feature is missing when using the 
maximally- localized Wannier basis. 

Table IVTll shows final occupations of the Wannier or- 
bitals varying with temperature and interaction strength. 
The values for (/3 oo,U = J=0) correspond to the LDA 
limit. Comparing with LDA values obtained in Ref. [93l |. 
the current LDA orbital polarization is not so severe in 
the Wannier description. This is due to the fact that in 
Ref. we used an empirical downfolding for the local 
DOS, merging most sulfur character with the Aig band. 
Since in the Wannier description the downfolding of sul- 
fur now also puts some weight on the Eg states, their 
filling is somehow increased. Nonetheless, turning on U 
does transfer sizable charge between the orbitals, hence 
the correlation effects envisioned in Ref. are indeed 
confirmed in the more elaborate Hamiltonian framework 
with WFs. However, due to the now resolved Aig-Egi 
hybridization, the charge transfer dominantly takes place 
between these two orbitals, leaving Eg2 as a mere "spec- 
tator". These interorbital charge transfers suggest that 
the FS of this material might actually be quite differ- 
ent than predicted by LDA, namely that the relative size 
of the various FS sheets may be significantly changed. 
A word of caution is in order however: the Luttinger 
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FIG. 17: (Color online) (a) Wannier DOS for BaVSs in the 
crystal-field basis, (b) corresponding k-integrated spectral 
function for BaVSa within LDA-hDMFT for 13=30 eV"^ 



theoremiS^ only constrains the total fc-space volume en- 
compassed by all sheets of the FS, and stipulates that 
it should correspond to one electron per vanadium for 
the present material, independently of interactions, in 
the metallic phase. There is no a priori theoretical re- 
lation between the volume of each individual sheet and 
the orbital populations as calculated above. Nonetheless, 
the reduction of the Aig population in favor of Eig may 
provide a hint that there is a corresponding shrinking of 
the Aig FS sheet. We now address in details whether 
this shrinking does occur and along which directions in 
reciprocal space. 

In order to investigate k- resolved effects of correlations, 
one must in principle determine the real-frequency self- 
energy. Since we are however mainly interested in the 
Fermi surface of the interacting system, we can extract 
the low-energy expansion of 'S{lu) from our QMC calcu- 
lation in the form: 

5ftS„„„, {u; + tO+) 2± 5RS„„„, (0) + (l - [Z-i]„„„,) lj + ■ ■ ■ , 

(28) 

where Z describes the matrix of QP weights. The QP 
dispersion relation is then obtained from the poles of the 
Green's function: 



det[wkl - Z (HKs(k) + 5ftS(0) - ^1)] = 



(29) 



Correspondingly, the FS in the interacting system is de- 
fined by 



det[At]l - HKs(k) - 5RS(0)] ^ 



(30) 



From Eqs. H29I30|I one understands that KS(0) pro- 
vides an energy shift to the LDA bands. The direction 
and magnitude of this shift depends at each fc-point on 
the amount of the contributing orbital character, since 
5ftS(0) is exphcitly orbital-dependent. Hence although 
our self-energy is explicitly k independent within single- 
site DMFT, one may still evaluate some k-depcndent ef- 
fects due to the explicit HKs(k) inclusion. 

From the converged self-energy matrix S(iw„) we de- 
rived 5RS(0) and Z via Pade approximation. Table IVlTll 
displays the two matrices. First, the symmetry of these 
matrices follows the earlier observations, hence there is 
some Aig-Egi coupling also reflected in the self-energy. 
The band renormalization cast into Z are roughly of the 
same order, with a slight maximum band narrowing close 
to a factor of two for the Egi band. Since we want to elu- 
cidate FS deformations due to correlations, looking for an 
explanation for the CDW instability, the KS(0) matrix 



TABLE VII: Band fillings for BaVSs from LDA+DMFT 
within the crystal-field Wannier basis. 



P (eV-^) 


U,J (eV) 


Aig 


Egl 


Eg2 


00 


0.0, 0.0 


0.59 


0.31 


0.10 


30 


0.0, 0.0 


0.58 


0.30 


0.12 


30 


3.5, 0.7 


0.41 


0.45 


0.14 
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TABLE VIII: KE(0) (in eV) and Z for BaVSa in the Wannier 
crystal-field basis from LDA+DMFT. 







KS(0) 






Z 






Alg 


Egl 


£^92 


Alg 


Egl 


Eg2 


Alg 


0.35 


-0.01 


0.00 


0.57 


-0.02 


0.00 


Egl 


-0.01 


-0.03 


0.00 


-0.02 


0.52 


0.00 


Eg2 


0.00 


0.00 


0.03 


0.00 


0.00 


0.58 



is of high relevance. From that, states with strong Aig 
character should considerably shifted upwards in energy, 
while for dominantly Eg states it depends on their sym- 
metry. Bands with Egi character should be shifted down, 
whereas those with Eg2 character should be shifted up. 
The latter discrepancy may be relevant to understand the 
opening of the gap at the MIT within the Eg states. Since 
the dispersion of the Aig-like band is highly anisotropic, 
i.e. ID-like, the rather strong shifting should result in 
a major FS deformation, invoking the possibility for an 
arising CDW instability. 

Fig. El displays the QP band structure according to 
Eq. H29|l close to the Fermi level. Using the hnearized 
self-energy should be valid within the small energy win- 
dow around s-p . Note first the overall narrowing of the 
bands since Z<1. The former statements concerning the 
respective shifts of the bands are accordingly reproduced. 
The electron pocket at the F point was identified as Eg2 
like (compare Fig. 112b ') and hence is now shifted up- 
wards. Roughly speaking, the system is getting rid of 
the Eg2 states when turning on correlations. Also the 
-Egl -like bands are considerably shifted downwards. In 
fact it is revealed in studies of the monoclinic phase be- 
low Tmit, that the internal Eg splitting seems indeed 
realized in the way outlined aboveiS^. Moreover, the off- 
diagonal Aig-Egi self-energy terms within LDA-I-DMFT 
lift the band degeneracy in the E-T-Z plane of the BZ. 
Originally it was arguediS^ within the LDA picture that 
the presence of the doubly-degenerated bands at this 
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FIG. 18: (Color online) k-resolved correlation effect for 
BaVSa. QP band structure from LDA-fDMFT (darkgreen 
dashed/dark dashed) in comparision with the Wannier band 
structure based on LDA (cyan/light gray) close to the Fermi 
level. 



zone boundary ensures the metallic properties of Cmc2i- 
BaVSa- Concerning the ^ig-like band along F-Z, there 
is a small shift to a lower kp, however not as strong 
as expected from Tab. I Villi This may be explained by 
the fact that the character of this band gains substan- 
tial Egl weight when approaching the Fermi level (to be 
observed in Fig. I15|l . Thus the states away from Sp are 
strongly shifted upwards but this shift weakens for the 
very low-energy regime. Recall that also in Ref. and 
recent angle-resolved photoemission (ARPES) measure- 
mentsiSSiiiS the shift to a lower kp{Aig) along F-Z was 
also not too strong. 

To proceed it is important to realize that studying only 
the single high-symmetry line along F-Z in order to re- 
veal a possible CDW instability may not be sufficient. In 
fact FS nesting can only be thoroughly investigated by 
taking into account the complete FS in the full BZ. With 
this aim we have plotted the LDA FS and the deformed 
LDA-fDMFT FS in Fig. [H In LDA the FS consists of 
two sheets which for the main part can truly be associ- 
ated to Alg and Eg. The ^ig-like one shows the expected 
strong ID behavior, i.e., corresponds to two main surfaces 
extending dominantly along {kx,ky), however with no real 
flattening. As already stated, the distance between these 
two surface parts is too large to account for the exper- 
imentally observed nesting. The second sheet incorpo- 
rates the Eg2 electron pocket at F as well as two Egi 
"pillars" on the h* axis (where also F-Y runs), extend- 
ing along c* and opening towards the zone boundaries. 
Within LDA+DMFT the Eg2 pocket shrinks and the Egi 
pillars thicken. Interestingly, the first sheet is indeed only 
little modified in the immediate neighborhood of the c* 
axis, however the further parts along a* tend to fiap down 
towards lower absolute values. Finally these latter parts 
are now not only within the distance of the experimental 
nesting vector qc=0.5c*, but are furthermore somehow 
more fiattened. That this shift indeed brings these outer 
parts very close to half filling may be seen in Fig. E51 
There the Wannier band-structure plot includes a line 
with non-high symmetry endpoints named "M/2" and 
"A/2". These latter points are halfway from F to M 
and halfway from Z to A (cf. Fig. I19() . A true quan- 
titative determination of the influence of the FS shape 
and of nesting properties on the CDW instability will 
require a full calculation of the Lindhard function us- 
ing the presently calculated self-energy. Although it is 
conceivable that a truly accurate account of the CDW 
transition for this material would require going beyond 
a k-indcpendcnt self-energy, we feel that the interorbital 
charge transfer and corresponding FS deformation found 
in the present work qualitatively points to the correct 
mechanism. 

Let us finally comment on the role of the Eg electrons 
in the MIT. It appears from our findings on the FS that 
the CDW in this material is in fact associated with the 
Alg electrons. The question then arises of the fate of the 
remaining Eg electrons through the MIT. Of course, from 
a many-body viewpoint, and in view of our present and 
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c. LDA+DMFT 






FIG. 19: (Color online) (a) 1. BZ of the Cmc2i structure, 
including the high-symmetry lines used to plot the band struc- 
tures. The tt* axis runs towards M, the b* axis towards Y. (b 
column) LDA Fermi surface, and (c column) QP Fermi sur- 
face from LDA-I-DMFT. The second row of the columns show 
the BZ along a* , the third row along b* . In the latter column 
the experimental CDW nesting vector qc is displayed. 



earlier resultsS^, one has to question the band-like char- 
acter of these states from early on. The very narrow Eg 
bandwidth causes a very low coherence scale for the cor- 
responding QPs. Thus above Tmit the Eg QPs have not 
yet reached their coherence scaleS^, which leaves them 
essentially localized already in the "metal". In a sense, 
this is a realisation of an "orbital selective" phasei^i 
at intermediate temperatures, i.e., a coexistence region 
of correlation-induced localized states and still itinerant 
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FIG. 20: (Color online) (a) LDA band structure (dark) and 
according Wannier bands (cyan/light gray) along a closed 
path in the 1. BZ including the line connecting "M/2" and 
"A/2" (see text), (b) pure (cyan/light gray) and renormal- 
ized (darkgreen dashed/dark dashed) Wannier bands. 



states for T>Tmit- In such a regime, the Aig electrons 
acquire a large scattering rate due to the presence of 
almost l ocali zed Eg degrees of freedom, as studied in 
Refs. |2nllTB . This is a likely explanation of the "bad- 
metal" behavior observed above the MITiSfliiiSIS^, and is 
also in good agreement with recent optical studiesii^. We 
also performed supporting two-site cluster-DMFT calcu- 
lations, in order to allow for explicit intersite V-V self- 
energy terms, revealing no essential differences to the out- 
lined picture. A closer investigation of the CDW state, 
i.e., the electronic structure of the system in the param- 
agnetic insulating regime below Tmit is however neces- 
sarji22iiflLii4, including their spin degree of freedom, to 
reveal more details about the role of the Eg elecrons. 



IV. CONCLUSION AND PROSPECTS 

In this paper, we have explored in detail the use of WFs 
as a flexible technique to perform electronic-structure 
calculations within the LDA-I-DMFT framework. WFs 
are useful for two different purposes. First, as a physi- 
cally meaningful way of defining the correlated orbitals 
to which a many-body treatment will be applied. Sec- 
ond, as a convenient choice of basis functions for inter- 
facing the many-body (DMFT) part of the calculation 
with virtually any kind of band-structure method. In this 
paper, three different methods have been used, namely 
a pseudopotential-based method (MBPP), an FLAPW 
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method and an LMTO-based method. We have ap- 
plied this approach to two transition-metal compounds. 
The first one, SrVOs, was chosen in view of its sim- 
phcity as a testing ground for performing a detailed 
comparison between two different Wannicr constructions: 
the maximally-localized- Wannicr-function construction 
(MLWF) and the N-th order muffin-tin-orbital (NMTO) 
method. For this simple material, very close agreement 
between the two methods was found. The second mate- 
rial, BaVSs, was chosen in view of its physical interest 
and of some open key questions. We have been able, 
in particular, to make new precise statements about the 
correlation-induced changes of the Fermi surface for this 
material, using our Wannicr based method. 

There are several directions in which the present work 
can be extended, and several open issues which need fur- 
ther attention. Here, we outline just a few of them. 
First, we emphasized above that the localisation prop- 
erties of the Wannicr functions depend of course on the 
energy range (or alternatively the set of bands) defining 
the subspace W in which the Wannier construction is 
performed. The correlated orbitals are then picked as a 
subset CCVV, in general smaller than W itself. In the 
actual calculations performed in this paper, the minimal 
choice C~W was made, associated for the two materials 
that were considered to exhibit low-energy bands with 
dominant t2g character. This of course, means that the 
correlated orbitals defining C had sizeable weights on lig- 
and atoms. In contrast, one may want to enlarge W (in- 
cluding in particular ligand bands) and define the corre- 
lated subset from Wannier orbitals which would be more 
localized on the transition-metal site. Exploring these 
various choices, and comparing them to other choices in 
which the correlated orbitals are not constructed from 
WFs (e.g are taken to be heads of LMTOs, or truncated 
WFs, or even atomic wavefunctions) is certainly worth 
further investigations, particularly in the context of late 
transition-metal oxides. An important related issue is 
the appropriate way of calculating the local-interaction 
(Hubbard) parameters from first principles, for each of 
the possible choices of correlated orbitals. Working with 
WFs will make it easier to address this issue in a manner 
which is independent of the underlying band-structure 
method and basis set used for performing the calcula- 
tions in practice. Finally, one should emphasize that the 
accuracy of the DMFT approximation depends on the 
choice made for the local orbitals defining C. In fact, 
it has been recently suggestedii^ that the local orbitals 
could be chosen in such a way as to make the DMFT 
approximation optimal (according to some criterion on 
the magnitude of the local interactions). 

Another line of development that we arc currently 
pursuing, is the practical implementation of the self- 
consistency over the charge density (and of total energy), 
along the lines of Appendix ^ within non-LMTO based 
electronic structure methods. A more remote perspec- 
tive for such developments would be the possibility of 
allowing for local structural relaxation for correlated ma- 



terials within LDA+DMFT. Additionally, the Wannier- 
based formalism is also well suited for the calculation of 
response functions (e.g. optics). We hope to be able to 
address these issues in future work. 



APPENDIX A: SELF-CONSISTENCY OVER THE 
CHARGE DENSITY IN LDA+DMFT 

In this appendix, we bricfiy discuss the implementa- 
tion of self-consistency over the charge density in the 
LDA-fDMFT framework. Until now, this has been im- 
plemented in practice only in the MTO^iiiLiiLiiS or 
KKRii^ context. Here, we discuss the use of other basis 
sets, with particularly the Wannier framework in mind. 

As outlined in section (|II A 4|) . charge self-consistency 
means that the KS cycle and the DMFT loop are ex- 
plicitly coupled, i.e., the charge density is calculated at 
the end of a DMFT cycle (including self-energy effects). 
From the new charge density, a new KS potential is ob- 
tained, and so on. 

The charge density is calculated from the full Green's 
function of the solid. Without yet introducing a specific 
basis set, it is given by Eq. ((T^ . namely: 




iujn + fJ-+ Vks - AS 



We shall find it convenient to split p(r) into 

p(r) = PKs(r) + [p{r) - PKs(r)] 
= PKs(r) + Ap(r) , 

with PKs{r) = (r Gks i". 



(Al) 

(A2) 
(A3) 



It is important to realize here that the demand for charge 
neutrality is not imposed on jOKs(r) but rather on p{r). 
The chemical potential /i must be therefore explicitly de- 
termined (at the end of a DMFT loop) in such a way 
that the total number of electrons is the correct one (i.e. 
charge neutrality is preserved): Ne=tT G=J drp{r). This 
value of the chemical potential will in general not be such 
that trG'KS= / c?rpKs(r) equals the total number of elec- 
trons Ne- This is quite natural, since the KS representa- 
tion of the charge density by independent KS wavefunc- 
tions no longer holds in the LDA-I-DMFT formalism. 

Formally, the charge-density correction Ap may be ex- 
pressed as 



Ap{r) 



r G - Gks r 



G 



KS 



Gj^g - G 



G 



Gks AEG 



(A4) 



In a concrete implementation, this equation must be 
written in the specific basis set of interest. We shall do 
this with the Wannier formalism in mind, i.e. use the 
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basis functions H19|l corresponding to the set W, which 
are unitarily related to the Bloch functions by 



\Wka 



E^-l^kJ- (A5) 

ue w 



Since in the Wannier basis set Gks and G are block diag- 
onal, they have common non-zero matrix elements only 
in the chosen subspace W. Let us define the operator 
AiV in W by its matrix elements: 

^^aJ' = l E GKS(k,iL^„)AI]^^,(zc^„)G™'a'(k,zc^„), 

n mm' 

(A6) 

Note that the indices aa' run over W, while the sum- 
mation over mm' runs only within the correlated subset 
C. Note also that the frequncy summation in this ex- 
pression converges rather quickly, since the terms in the 
sum decay as 1 /w^ at large frequencies (hence no conver- 
gence factor e*'^""^ is needed in the sum). The desired 
charge-density correction therefore reads: 

Ap(r) = E ('■l^ko) Aivfj, {w^^, |r) (A7) 

kaa' 

Alternatively, IjATp can also be written in the Bloch basis 
as: 

Mr)-E E V'k.W^L'WATV^::?, (A8) 
with ANi':;}^Y.Ul^^^N^3u!^l' (A9) 



The KS part of the charge density is easily calculated 
within the given band-structure code: 



Pks 



(r)= E IVv(r) 



(AlO) 



Note the difference between HA8|I and HA10|I . as in (|A8p 
additional terms, off-diagonal in the band indices, con- 
tribute to Ap(r). In this respect it proves to be con- 
venient to introduce the density matrix of KS orbitals: 



^^'^(r) = V'k.(rW'L'(r), 



(All) 



The two contributions to the charge-density may now be 
compactly written as: 



Pks 



(r) = E (^) ' 



kiy 



Ap(r) = E*"^-^ 



(A12) 



(A13) 



Hence finally, we can write the full charge density, i.e., 
eq. (|A1|) . in the following explicit form: 

pir) ^ E D';!:l{r)ANl^} + E ©(m " e^.)Dl}$ir) . 

(A14) 



Generalization to finite-T DFT is straightforward. Real- 
ize that in the Wannier implementation, the first (dou- 
ble) sum in ljA14p would only run over Bloch bands in 
the W set (since only there AA^ is non-zero), while the 
second sum runs over all filled Bloch states. Again, it is 
important that the chemical potential has been correctly 
updated at the end of the DMFT loop, so that the full 
charge density (|A14(I sums up to the correct total number 
of electrons, hence insuring charge neutrality. 

Expression l|A14(l has in fact a general degree of valid- 
ity, not limited to the Wannier implementation, provided 
the matrix AA^, originally calculated in the basis set in 
which the DMFT calculation has been performed, is cor- 
rectly transformed to the Bloch basis set. 



APPENDIX B: FREE-ENERGY FUNCTIONAL 
AND TOTAL ENERGY 

For the sake of completeness, we briefly summarize in 
this appendix how the LDA-I-DMFT formalism for elec- 
tronic structure, described in Sec.^can be derived from 
a free-energy functional. As pointed out by Savrasov and 
Kotlias^, a ( "spectral-density-" ) functional of both the 
electron charge density p(v) and the on-site Green's func- 
tion in the correlated subset Gj^^, , can be constructed 
for this purpose. Let us emphasize that these are inde- 
pendent quantities: because Gioc is restricted to on-site 
components and to the correlated subset C, the charge 
density p(y) cannot be reconstructed from it. A com- 
pact formula for the total energy can also be obtained 
from this functional formulation£i^2iiSi, given at the end 
of this appendix. 

The functional is constructed by introducing 
source terms T4cs(r)-^K:xt(r) and AEmm'(itJ„), 
coupHng to the operators ^'^{v)i\;{v) and 
Et X™(r-R-T)^(r, t)^\v\ r')x™'(r'-R-T), re- 
spectively. It reads: 

0[p,G5r™';"^"KS,AI]™„,] = (Bl) 
1 



V2 

--tr In icjn + + — - Vks - AS 
- y" rfr (yKs(r) - T4xt(r)) p(r) - tr[GiocAS] 



T 



r - r' 



($imp[G 



*Dc[G5^m']) 



(B2) 



In this expression, $iinp is the Luttinger-Ward functional 
of the quantum impurity model, and ^dc a correspond- 
ing functional generating the double-counting correction. 
Note that, as usual, it is understood that AE in the log- 
arithm term of (|B2p represents the self-energy correction 
on the lattice. 

Variations of this functional with respect to the sources 
(517/(5 T4cs=0 and (5ri/(5Emm'=0 yield the standard expres- 
sion of the charge density and local Green's function in 
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terms of the full Green's function in the solid, Eq. ||2Jl. 
The Legendre multiplier functions Vks and AE can be 
eliminated in terms of p and Gioc, so that a functional 
only of the local observables is obtained: 

r[p, Gioc] = n [p, G'Za';VKs[p, Giod, As[p, dod] . 

(B3) 

Extremalization of this functional with respect to p 
{6r/Sp=0) and Gioc ((^r/JGioc—O) yields the expres- 
sion of the KS potential and self-energy correction at 
self-consistency, i.e., the self-consistency conditions (over 
both the local projected Green's function and the charge 
density) of the LDA-fDMFT formalism. 

Using the above expressions, the free-energy can be 
written as: 

n = fluFT + ti' In Gj^g — tr In G^^ — tr [GiocSjmp] 

+tr [G IocSdc] + J2 (^i-P - *Dc) ■ (B4) 



In this expression, J^dft is the usual density-functional 
theory expression of the free-energy (calculated at the 
self-consistent LDA+DMFT charge density, however). 
Taking the zero-temperature limit of this expression leads 
to the following exprcssion^iiSS of the total energy at 
T=0: 



-Elda+dmft = Edpt + {Hu) — Euc 



kact' 



/KS 



(B5) 



Note that a, a' cover the full electronic Hilbert space. 
Finally, it should be understood that in deriving the 
LDA-fDMFT equations from this free-energy function- 
als, we have implicitly assumed that the orbitals Xm 
defining the correlated subset C were kept fixed. In prac- 
tice however, one may want also to optimze these or- 
bitals, e.g. by minimizing the free-energy Q. Further- 
more, we emphasize that in the implementation of the 
charge self-consistency described in Appendix 1X1 it was 
explicitly assumed that these orbitals were recalculated 
as Wannier functions associated with the new KS poten- 
tial. This is done at each stage of the iteration over the 
charge density. If instead the local orbitals are frozen 
(e.g. from the LDA potential), then they are no longer 
unitarily related to the new set of KS orbitals correspond- 
ing to the new potential, and the formulas derived in 
Appendix 1X1 have to be appropriately reconsidered. 



constructed for Sr((4s4p4(i)/Ba(5s5p5(i), V{3s3p3d) and 
0((2s2p2(i)/S((3s3p3(i). The partial-core correction was 
used in all constructions. In the crystal calculations, lo- 
calized functions were introduced for all valence states. 
Thereby, the localized functions are atomic pseudo wave- 
functions, either multiplied with a cut-off function (sd) or 
minus a spherical Bessel function (p) . The cut-off radius 
for the V{3d) local functions was chosen as 2.0 a.u.. The 
plane-wave cutoff energy E^^^ was 24 Ryd for SrVOs 
and 20 Ryd for BaVSa, and the Perdew-Wang exchange- 
correlation functionaliSi was used in all calculations. 

In the FLAPW calculation for SrVOs, the atomic- 
sphere radii for Sr, V, and O were chosen to be 2.5, 2.0, 
and 1.55 a.u., respectively and the FLAPW basis size 
was set to include all plane waves up to £'p™'''=20.25 
Ryd . The local orbitals sector of the basis was used 
to include the high lying V(3s), V(3p), Sr(4s), Sr(4p) 
and 0(2s) core states as valence states in the calcula- 
tions. Additional V{3d) local orbitals were used to relax 
the linearization of the transition metal d-bands. The 
Hedin-Lundquist exchange-correlation functionalist was 
used in the calculation. 

For the self-consistent LMTO calculation of the effec- 
tive LDA potential for SrVOa we have used the prescrip- 
tion described in Ref. . The radii of potential spheres 
were chosen to be 3.46 a.u., 2.3 a.u. and 2.04 a.u. for 
Sr, V and O atoms, respectively. Twelve empty spheres 
were introduced above the octahedron edges with the 
radii 1.01 a.u. The LMTO basis set used in the cal- 
culation were s{p)d{f), spd, {s)p{d) and s{p){d) for Sr, 
V, O and empty spheres, respectively, (l) means that the 
Z-partial waves were downfolded within TB-LMTO-ASA. 
The von Barth-Hedin exchange-correlation functionalist 
was used in the calculation, which is very similar to the 
Hedin-Lundquist parametrization. These elder i?xc rep- 
resentations utilize the same interpolation formula. 

The construction of the MLWFs from the disentangling 
procedure for BaVSa was described in section UlIB 31 In 
the case of SrVOs, no energy windows have to be de- 
fined, since the t2g-likc bands form an isolated group of 
bands. The NMTO-WFs for SrVOs were obtained by 
imposing that for a <2g orbital, the <2g character on any 
other V site must vanish. We used N=2 and chose the 
energies for which this quadratic MTO set is complete 
as: eo=:-1.29 eV, ei=0 eV and £2 = 1-61 eV with respect 
to Fermi level. Finally, the NMTO set was symmetrically 
orthonormalized. 



APPENDIX C: TECHNICAL DETAILS 

In the following we provide further details on the band- 
structure calculations performed in this work. Con- 
cerning the computations using the MBPP code, due 
to the significant ionic character of the treated com- 
pounds, the semicore states of Ba, Sr and V were 
treated as valence. Hence the pseudopotentials were 
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